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Short path length pQCD corrections to energy 
loss in the quark gluon plasma. 

by Isobel Kolbe 

Recent surprising discoveries of collective behaviour of low-pT particles in pA collisions 
at LHC hint at the creation of a hot, fluid-like QGP medium. The seemingly conflicting 
measurements of non-zero particle correlations and RpA that appears to be consistent 
with unity demand a more careful analysis of the mechanisms at work in such ostensibly 
minuscule systems. We study the way in which energy is dissipated in the QGP created 
in pA collisions by calculating, in pQCD, the short separation distance corrections to the 
well-known DGLV energy loss formulae that have produced excellent predictions for A A 
collisions. We find that, shockingly, due to the large formation time (compared to the 
1/p Debye screening length) assumption that was used in the original DGLV calculation, 
a highly non-trivial cancellation of correction terms results in a null short path length 
correction to the DGLV energy loss formula. We investigate the effect of relaxing the large 
formation time assumption in the final stages of the calculation - doing so throughout 
the calculation adds immense calculational complexity - and find, since the separation 
distance between production and scattering centre is integrated over from 0 to oo, > 100% 
corrections, even in the large path length approximation employed by DGLV. 
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Chapter 1 


Introduction 


1.1 Motivation 

While the SKA will soon be pushing the boundaries of South African physics, particle 
physics rivals cosmology on the global stage as one of the most fruitful field of discovery 
physics, and it is truly the heavy ion particle physicists that are pulling down the veil of 
obscurity surrounding the underlying structure of the universe humanity finds itself in. 
The heavy ion community has all but universally accepted that a new state of matter is 
routinely created in colossal particle colliders such as the Large Hadron Collider (LHC) 
at CERN in Geneva [1], Switzerland and the Relativistic Heavy Ion Collider (RHIC) [2] 
at Brookhaven National Labratory (BNL) in the United States. This new state of matter, 
called the Quark Gluon Plasma (QGP), is thought to be a deconfined state of the most 
fundamental constituents of matter, quarks and gluons, that is brought about by the 
extreme temperatures and energy densities created by colliding heavy nuclei like gold at 
RHIC and lead at LHC at ultra-relativistic energies. The QGP created in accelerators 
offers unique insight into the nature of the basic components of matter and it is therefore 
imperative that mankind investigate its properties. To date it is known that the QGP is 
extremely short lived, with a hadronization time on the order of 10 fm/c [1], an attribute 
that excludes the use of any externally produced probes to discover the properties of this 
medium. In order to characterize the QGP then, one must rely on self-generated probes 
- probes that are produced along with the medium. 

In a heavy ion collision (hereafter ‘AA collisions’), thousands of particles are produced 
(~ 7000 in central AuAu at RHIC [3], ~ 10 000 in central PbPb at the LHC [4]). The 
legion of particles are subjected to an intricate process that starts with a difficult-to- 
determine initial nuclear parton distribution, followed by ‘melting’ into the QGP which 
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itself evolves obscurely before hadronizing into particles that finally move out to a de¬ 
tector. Particle physicists are limited to data taken meters away from a medium that 
is only femto-meters across and must therefore rely on evidence from a conglomerate of 
phenomena in order to make scientifically sound deductions regarding the nature of the 
QGP phase of the collision. In essence, what can be learnt from a detector is only the 
type and momentum of the particles produced. Theorists must then attempt to inter¬ 
pret the measured distributions of the particles and their momenta. We will focus in the 
Introduction on two major signatures of the QGP that appear in the momentum distri¬ 
butions of the particles. The first is known as ‘flow’, which indicates collective behaviour, 
while the second is known as ‘quenching’ or ‘particle suppression’, which indicates that 
particles have lost energy as they propagated through the medium. The dual phenomena 
of elliptic (and higher order) flow of soft (or low transverse momentum, pj') observables 
and the quenching of hard (or high pt) observables are expected wherever a hot medium 
is created [1]. 

An early solution to the problem of controlling the initial state was to simply attempt to 
measure it; in principle one should be able to determine the properties of the initial state 
if one can probe the nucleus without creating a medium. Since it was inconceivable that 
hot nuclear matter could be created in a collision between a heavy ion and a proton, RHIG 
performed a series of experiments involving different combinations of protons, deuterons 
and a range of heavier nuclei [5] from which it was concluded that RpA ~ 1 (or RuAu), sug¬ 
gesting that there were no hot nuclear matter effects present in these mixed experiments 
and that the results served as a decent control measure with which to carefully analyse 
AA data. Following suit, the LHG then performed a p-Pb (proton-lead) run in order to 
reproduce these results, but, shockingly, discovered signatures of hot QGD matter present 
in their pA experiments. Most remarkably, eight particle correlations seen by GMS [6] 
(shown in Figure 1.1.1), but also measurements of U 2 , and even (the Fourier modes 
of particle correlations) at ATLAS [7]. 

The presence of what appears to be collective behaviour in seemingly impossibly small 
systems (as any system of interacting particles was expected to be in a pPb or dAu 
experiment) has raised a number of uncomfortable questions: Firstly, collective behaviour 
is well understood within the framework of viscous hydrodynamics [3] , but hydrodynamics 
is a statistical approach and effectively assumes an infinitely large system of particles in 
order to describe their collective behaviour [8] ; Is a system of no more than ~ 2 fm large 
enough for the assumptions in hydrodynamics to hold? Secondly, could it be possible 
that the number of interacting nucleons does not always scale in the same way? That 
is, is it clear that the scaling is always by the number of participants (the sum of the 
particles entering into the collision) or the number of binary collisions (each participant 
might collide more than once)? So do we have legitimate confidence in constraints placed 
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Figure 1.1.1: Multi-particle correlations in pPb experiments at the CMS detector [6]. 

on the extent of the effects of geometry and nucleon distributions from In light 

of the fact that centrality and multiplicity in pPb show a much broader correlation (left 
hand panel of Figure 1.2.2), can we even say that we have determined the centrality, 
and with it the number of participants, in a pPb collision clearly? In short this second 
question boils down to asking whether or not measuring RpA is as simple as originally 
thought, whether the measurement can in fact be made reliably, and if not, the extent 
to which the Raa results are under experimental control. A third problem, though not 
new and certainly not unique to the pA case: Virtually all QGP phenomena can be 
described by either low pT, strong coupling physics like hydrodynamics [ 8 ], or high px, 
weak coupling physics like perturbative QCD (pQCD)[9], but if we believe that both 
formalisms describe nature accurately, then the coupling must run, and yet it is still by 
no means clear how the coupling is to run. We must attempt to discover whether or not 
high momentum observables can be understood in pA experiments within the framework 
of perturbative QCD (pQCD), as they have been in AA collisions [10]. If a QGP does 
exist, the experiments must show energy loss. 

While not monopolizing QGP observables, V2 and Raa offer important insight into col¬ 
lective behaviour and jet quenching, the careful study of which is the key to gleaning 
information about the QGP. However, the experimental data are mute without a rigor¬ 
ous theoretical understanding of the underlying mechanisms that produce the observed 
phenomena. Jet quenching has been exceptionally well described by a number of energy 
loss formalisms, but in this work the GLV (Gyulassy, Levai, Vitev) setup [9, 11, 12] was 
followed as it is pedagogically insightful. The GLV description of radiative energy loss 
has been immensely successful in describing the energy loss of partons moving through a 
QGP that is created in very central A A collisions, as can be seen in Figure 1 . 1 . 2 . Figure 
1 . 1.2 shows that initial state effects are well under control since we see no suppression 
of the (colorless) photons. Furthermore, particles with very different masses, the and 
r] (~ 135 MeV and ~ 548 MeV, respectively), show similar suppression, indicating that 
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Figure 1.1.2; 


Raa for AuAu at RHIC compared to GLV energy loss for 7, 7r° and ij 
[13, 14]. 


energy is lost partonically rather than hadronically - that is, the constituent quarks have 
already experienced energy loss by the time they are able to hadronize into observable 
hadrons. 

The GLV predictions are the result of a number of calculations and corrections to previous 
works [9, 11, 12], resulting in a sophisticated calculation that takes into account finite time 
effects (looking at the energy loss of a particle created at t = 0 inside the medium rather 
than at t = — 00 ) and ultimately also generalizing to massive quarks, the DGLV (Djordje- 
vic GLV) result [10]. The GLV result makes use of the Gyulassy-Wang model [15], which 
describes in-medium interactions by considering a number of static scattering centres and 
then looking at a very high momentum parton (a quark or a gluon) that scatters off 
the centres and is thereby engendered to radiate a gluon. However, the calculations rely 
heavily on a simplification that arises when one considers only large separation distances 
between the static scattering centres. Gonsidering only large separation distances means 
that the calculation is only valid for large system sizes (as the separation distance must 
necessarily be smaller than the system size in order for any scattering to occur). There¬ 
fore, even though the result obtained by DGLV [10] is fully general for massive quarks 
and indeed depends linearly on the size of the system, it is in fact not valid for small 
system sizes since the derivation of the result makes repeated use of a large separation 
distance assumption. The monumental success of energy loss studies in AA collisions, 
that have been extended even beyond the DGLV result to include dynamical scattering 
centres [16, 17], provide an excellent starting point for carrying out a short path length 
extension. 

A need has therefore arisen for a quantitative analysis of energy loss in pA collisions — 
not only to determine whether or not a droplet of QGP exists in pA experiments, but 
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also to ensure that all relevant effects are properly under control. In order then to allow 
the field of high energy particle physics — in particular heavy ion physics, to advance, it 
is vital that the effects seen in pA be understood wholly. A precise statement regarding 
whether the coupling of the medium is strong or weak, the applicability of competing 
theories and the extent to which current interpretations of RpA ~ 1 hold, cannot be made 
without a falsifiable energy loss calculation. It is exactly the absence of such a prediction 
that we aim to address. 
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Figure 1.2.1; Impact parameter and centrality, (a) The reaction plane defined by 
HP and impact parameter b.[3]. (b)Yield as a function of total charge collected in the 

beam-beam counters [18]. 

1.2 Relevant Concepts and Current Problems 

In this section I review the various concepts and current difficulties that inspired the 
execution of the calculation in this dissertation. Although some of the topics below do 
not enter explicitly into the calculation, they serve as important aspects of the context in 
which the main calculation is performed. 

Centrality and Mnltiplicity 

Figure 1.2.1a shows the definition of the impact parameter b which is the spatial vector 
between the centres of the two colliding nuclei. The magnitude of the impact parameter 
defines the ‘centrality’ of an event while its direction defines the event raction plane. For 
historical reasons, centrality is given as a percentage in such a way that the most central 
events (smallest impact parameter) have the lowest percentage while more peripheral 
events (large impact parameter) have a larger percentage. Figure 1.2.1b is a plot of the 
probability of an event occurring in a given centrality bin, with most of the events that 
occur in a collider being peripheral (95% centrality) and only very few events in the 0 — 5% 
centrality bin. It is impossible to know a priori what the centrality of any given event is 
and centrality must therefore be determined on an ‘event-by-event’ basis. To determine 
centrality, experimentalists often use an observable called ‘multiplicity’^ which is simply 
the number of charged hadrons (charged mesons or baryons) measured in the detectors. 
In A A collisions, the centrality is directly related to the multiplicity, as shown in the right 
hand panel of Figure 1.2.2, which is a plot of the relationship between multiplicity and 
the number of participants as calculated in a Glauber Monte-Carlo simulation [4]. 

^Other determinations of centrality come from, e.g. Zero Degree Calorimeters. 
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Figure 1.2.2: The relationship between the number of participating nucleons and the 
multiplicity in pA collisions (on the left) and AA collisions (on the right). [4] 


The number of participant nucleons is closely related to the centrality (and the impact 
parameter), something that is easily seen geometrically. In the right hand panel of Figure 
1.2.2, we clearly see that in AA, the centrality is very closely related to the multiplicity - 
an extremely useful fact since centrality is an important quantity theoretically but must 
be measured indirectly experimentally through measurements of multiplicity. However, 
in the left hand panel, we see a broadening of the correlation, making it much more 
difficult to accurately relate an experimental measurement of multiplicity in pA to impact 
parameter. 

The poor correspondence between centrality and multiplicity in pA is often bypassed by 
considering ‘minimum bias’ results in which experimental observables are not measured 
for a given centrality. Instead, a kind of average over centrality bins is taken of all events, 
giving a result that is not centrality dependent. This smearing of centrality dependence is 
problematic in searches for the creation of ‘droplets’ of QGP in pA experiments because 
it is only expected that the QGP be created in very central pA collisions - since it is only 
in these events that one can expect a high enough energy density and enough particles to 
form an interacting system [8]. However, the converse solution is equally ambiguous; if one 
were to consider only high multiplicity pA events, the small number of particles involved 
introduces a selection bias that enters into the uncertainty related to a high pT jet - one 
cannot say whether a large number of low momentum particles in a given region are the 
decay products of one very high momentum particle (a jet) or are simply a collection of 
soft particles collected in one place due to some geometrical factor. 

We must therefore be clear when discussing multiplicity and centrality: In AA exper¬ 
iments, the terms are interchangeable, but they must be treated with more care when 
discussing pA experiments. 
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Figure 1.2.3: A minimum bias Rppb result (left hand panel) vs. the same evaluation 
for different centralities (right hand panel) at forward rapidity presented by the ATLAS 
collaboration at Hard Probes 2015 [19] 


Raa and R^a 


Raa is an experimental observable that compares the particle yields in an AA collision 
to those in pp collisions and can be plotted as a function of, for example, momentum or 
rapidity. Equation (1.2.1) shows how RaA) as a function of transverse momentum pp, is 
expressed in terms of the differential particle number seen in AA collisions {dN aa/ dpx), 
the differential cross section dapp/dpT in pp collisions and the average thickness function^ 

{TAA)[m- 


Raa{pt) 


1 dNAA/dpT 
{Taa) dapp/dpT 


( 1 . 2 . 1 ) 


An Raa = 1 should indicate then that not only is the initial state well understood, but 
also that there are no final state (medium induced) effects, while an Raa 7 ^ 1 should in 
principle speak to the presence of a medium, provided that all unrelated effects are well 
under control. One such effect that is not under control is the non-trivial distribution of 
nucleons within a nucleus which directly affects the counting of {Taa) while another is 
the afore-mentioned determination of multiplicity. The ambiguities associated with RpA 
are illustrated in Figure 1.2.3 which was presented by the ATLAS collaboration [4] and 
clearly shows the difficulties in interpretation connected to a claim that RpA ~ 1- 

The interpretational complexities associated with any RpA results require assurances that 
energy loss in a small systems is well under control - if one wants to claim that RpA ~ 1 
and that the scaling and initial state effects are therefore well under control, one must be 
absolutely sure that the effects governing energy loss are well understood. Furthermore, 


^Integrating over the average thickness function gives the average number of binary collisions. 
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the presence of multi-particle correlations (Figure 1.1.1) [6], although by no means suffi¬ 
cient support for the existence of hot QCD matter, does prompt the search for a QGP in 
these small systems. 

All successful energy loss calculations result in a (positive power) dependence on the size 
of the system and the expectation must then be that RpA measurements will yield small 
results even in the presence of a hot QCD medium. In light of the additional experimental 
difficulties, it is crucial that energy loss for short path lengths be understood from a 
theoretical standpoint in order to properly interpret RpA measurements. 

Other Energy-loss formalisms 

The current work is performed within the GLV framework, but I will discuss briefly here 
the features of some of the other formalisms that have been developed. 

There are four major formalisms that have been developed for computing partonic energy 
loss in a hot QCD medium using perturbative QCD (pQCD), each attempting to solve the 
problem using a different set of assumptions. The first radiative energy loss calculations 
were performed by BDMPS (Baier, Dokshitzer, Mueller, Peigne and Schiff) [21-24] in 
which multiple soft gluon scattering was considered. The major break from multiple soft 
gluon scattering was made by GLV in a series of papers [9, 11, 15] in which they developed 
a formalism to handle hard scattering in a thin plasma (where only a few scatterings are 
expected to occur). The ‘Reaction Operator’ formalism that GLV developed led to two 
more methods; ASW (Armesto, Salgado and Wiedemann) [25] and the ’Higher Twist’ [26] 
formalism. An alternative method was developed by AMY (Arnold, Moore, Yaffe) which 
is an effective Hard Thermal Loop theory [27]. Djordjevic and Heinz [17] made the next 
major leap in energy loss calculations by considering dynamical scattering centres [16] and 
it has recently been shown [28] that indeed, a dynamical approach is necessary to describe 
data more accurately. However, the static case provided theorists with excellent intuition 
in AA studies and therefore constitutes a useful first step that may provide insight with 
which to proceed. 

The analysis of energy loss in AA collisions is exhaustive and extensive and has in recent 
years been developed to an extremely sophisticated level [28, 29] taking into account a 
range of different effects including finite time considerations, both collisional and radiative 
energy loss, energy loss in both static and dynamic media, finite magnetic mass effects and 
the running of the coupling. In the near future it will be the task of theorists to extend 
this range of studies to short path lengths. We endeavour to lead this new analysis by 
starting, as was done in AA, with the static scattering case and generalizing the DGLV 
result for heavy quarks to include short separation distances. 



Chapter 2 


Formalism and Setup 


2.1 Notation and Conventions 


In this dissertation, for consistency with both [10, 12], I have used the following notation 
for vectors: 

• p: Transverse 2D vectors. 

• p = {pz,p)'- 3D vectors 

• p = {p^,p) = Ip^ + p^,p^ ~ P^jP]- four vectors in Minkowski and Light Cone 
coordinates respectively. 

The calculations are performed in Minkowski space, but we will use many values derived 
in light cone coordinates. The dot products are, in ‘mostly-minus’ Minkowski and light 
cone respectively. 


p.k = p^k^ - p^k^ - p • k 



P • k, 


( 2 . 1 . 1 ) 


and the transformation rule 


= hiP'^ +P~ 

= 5(p+ -p- 


p-< 


_ pU _j_ 


= P 


_ 0 z 

p = p — p 


p = p- 


Since the calculation is done in QCD, there will necessarily be color factors. Notations 
differ and so, in an attempt to be absolutely clear, the color exchanges are handled using 
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M 


1,0,0 




k,c 


/' 
1^0 / 


Z I 




P 


Figure 2.1.1: ^ i ^ o , o .[ 15 ] 


the applicable SU{Nc) generator Ta{n) in the dn dimensional representation of the target. 
That is dn refers either to the dimension of the fundamental representation {d{N) in [30]) 
in which the quarks live, or the dimension of the adjoint representation (d(G) in [30]) in 
which the gluons live. The generators Ta{n) are traceless, but obey ([10, 30]) 

Tt {Ta{i)n{j)) = 5,j6abC2{i)^ ^ TV(t“t^) = 5ab^,C2{r), (2.1.2) 

dA d[G) 

from the definition of the Casimir operator (in DGLV notation on the left and Peskin 
notation on the right). Therefore, in the notation used in this calculation (we will follow 
DGLV), dA is the dimension of the adjoint representation and C 2 {i) the Gasimir operator 
of whichever representation (adjoint or fundamental) has dimension di. In fact, when 
the color algebra is done explicitly, we will use the simplifying notation a = ta for the 
generator of the representation of the parton, with aa = CjiIr. Although these numbers 
are known (and can be found in [30]), the constants will be kept in algebraic form for 
clarity. 

A short hand for energy ratios will prove useful notationally as it has in the DGLV 
calculation. Gonsider therefore the diagram 2.1.1 defining the relevant momenta. The 
energy of the radiated gluon (with transverse momentum k) is given by w ~ = 

xE = xp^, with X the fraction of energy (or momentum) carried away by the radiated 
gluon and E the energy of the hard parton. We will then use 


xE~^ xP+ 


U) 


2 

k 2 


UlQ = 


Ui = 




2u} 

( 1 <^ - 
2oj 

(k - Ti - q, 


2lo 


mi + 


— 


2 cj 


(2.1.3) 









Formalism 


12 


Lastly, we will use a number of relations in the color manipulations. These are from [30] 


[a, b] = ab — ba = 
tr[a\ = 0 

facdfbcd ^ C2{G)6‘^^ 
aa = C 2 (r)l 
tr[ab] = C(r)5“^ 


(2.1.4) 


1 


r^^bc = -iC2{G)a 

rj. ,rr rr \ G2{T)dT 

TliJa^Ta^) = -- 6, 


ai,a2 
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2.2 Assumptions 

The present calculation is performed within the GLV formalism, the calculational details 
of which are discussed in Section 2.3. Within this formalism, a number of assumptions 
are made that affect both the ease of the calculation and the physical interpretation of 
the result. In this section I discuss these assumptions briefly. 


The Eikonal Approximation (High energy) 


In the eikonal approximation, the largest scale of the problem is the original energy of the 
hard parton. The momentum vectors that are relevant here are derived in section A.3 to 
be 


k = 

P = 


ml + k 

’ xP+ 


(1 - x)P^ 


m 2 +k2 
’ {l-x)P+ 



q = [q'^,q ,q]- 


( 2 . 2 . 1 ) 


Now, E~^ the largest scale of the problem means that P'’“ S> P . The eikonal approxima¬ 
tion therefore leads to the relation [12] 

= Wo ~ (2.2.2) 


Soft gluon and soft rescattering 

The radiated gluon carries away a momentum fraction x of the momentum of the parent 
parton. The soft gluon approximation, x <C 1, allows for a number of simplifications. 
(Note that, although x is a small number, E~^ is so large that we may assume xE~^ ^ jkj.) 

First, the eikonal approximation allows us to assume that the source current that describes 
the hard production of the parent parton varies slowly with momentum p. The scale 
on which this variation occurs is the Debye screened length the range of which is 
determined by the momentum transfers from the scattering potential. Therefore we have 
that, in light of both the soft radiation and eikonal approximations. 
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We can also say the following because of x <C 1: 

e,{k){p-q + kr^^ 

k-{p-q)^k-{p)^ — (2.2.4) 

We may also now extend equation (2.2.2) to read 

£1+ > = Wo ~ ■ (2.2.5) 


Impact parameter 


Here we assume that the impact parameter (describing the overlap of the nuclei in the 
collision) varies over a large transverse area A± relative to the interaction area 1 ///^. This 
assumption allows for a simplification to occur when performing the ensemble average 
over initial states. The impact parameter average reduces to 




which allows the ensemble average over the phase factor to become 

^2 




( 2 . 2 . 6 ) 


(2.2.7) 


It is worth noting that the ensemble average over the phase factor containing the impact 
parameter can still be reduced in this manner in pA since the transverse area over which 
the impact parameter varies is still large in comparison to the reaction area: the reaction 
area is ~ 1/(0.5GeH)^ ~ 0.1 fm^ while the overlap is, in central pPb which we 

consider, 7r(l/m)^ ~ 3fm^ 


Large formation time Ui -C /ii 


A finite formation time [15], 


T{k) ~ 


h 

AE{k) 


rsj 


2w 2 


( 2 . 2 . 8 ) 


for 9 = k/w, the angle at which the gluon is radiated, leads to a destructive interference 
phenomenon for particles emitted at small angles (with large momenta). The effect is 
known as the Landau-Pomeranchuck-Migdal (LPM) effect [31, 32], or coherent radiation 
limit, and causes a strong suppression of induced gluon bremsstrahlung [11] (not including 
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the ‘self-quenching’ phenomenon discussed in section 2.3) because a particle will scatter 
multiple times before radiating. LPM suppression is reduced in the case of a massive 
quark [10] because the formation time of a gluon emitted from a massive quark is reduced 
due to the non-zero mass effects. 


As in the calculation performed by [10], we will also assume that the formation time is 
much larger than the Debye screened length. 


/r 3> Wj 


2uj 


1 

r 


(2.2.9) 


This assumption only featured lightly in the computation of the DGLV result, but it will 
play a crucial simplifying role in the calculation presented in this dissertation and I will 
therefore elaborate on it briefly. 

By assuming that the formation time is large compared to the inverse Debye screened 
length, one is simply restating the GW model [15] - partons scatter off Debye screened 
scattering centres and therefore cannot resolve any structure of the scattering centre that 
might affect the interaction in a complicated manner. One expects to find that the energy 
loss due to radiation goes to zero for separation distances (separation between scattering 
centres) going to zero. However, if the formation length is larger than the mean free path 
one enters the LPM regime where a particle must scatter multiple times before radiating. 
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Mj 



Figure 2.3.1: and (g) .^o-[15] 


2.3 Formalism 

The calculation presented in the subsequent sections is a short path-length generalization 
of the result obtained by [10]. Therefore, the entire calculation is performed within the 
GLV formalism, based on the Gyulassy and Wang model [15] , then developed in a number 
of papers in which various assumptions were refined and corrected [9, 11, 12], and finally 
generalized by Magdalena Djordjevic in a paper [10] in which the GLV energy loss formula 
was generalized for heavy quarks and whose results are reproduced and added to in this 
dissertation. In this section I briefly review the GW (Gyulassy-Wang) model, following 
[9, 11, 12, 15], 

The GLV model considers the induced gluon radiation of a hard parton as it interacts 
with a random color field produced by a color neutral ensemble of static partons [15]. 
The setup therefore requires a means to handle (1) the production of a highly off-shell 
hard parton, (2) the vacuum radiation associated with such a hard parton and (3) the 
interactions with the medium that stimulate the emission of a soft gluon. Since we work 
in the high energy limit with a low momentum transfer, the spin characteristics of the 
parent parton can be ignored and we therefore deal with scalar QGD, a theory that is 
similar to scaler Quantum Electrodynamics (sQED), with the difference residing only in 
the addition of a treatment of color exchange. The setup below is therefore analogous to 
the sQED setup treated in [33] . 

Eor (1), consider the production (at finite time, in contrast to a jet that has been prepared 
in the infinite past as in [34]) of a hard, highly off-shell parton with momentum p. The 
parton is created along with the medium and is therefore (initially) localized at Xq = 
(0,x)^. The amplitude is assumed to vary slowly with p and can be written as [12] 


= ij{p)FP^nR 


( 2 . 3 . 1 ) 
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Even in the absence of a medium, a highly off-shell parton will radiate gluons, thereby 
losing energy (‘self-quenching’) and softening the spectrum of hard jets. To first order in 
gs, the softening of the spectrum of hard jets implied by equation (2.3.1) is effected as 
illustrated in Figure 2.3.1, and given (in the eikonal limit), by [10, 15] 


^vac - ( 2igs) ^2 + ^2 + 

We can therefore express the amplitude for a jet and soft gluon in the final state as [10] 

= -^J ^vac = iJ{p + k)F^^^^'>^°{iga){2p + k)^e^{k)i/S.M{p + k)c 

« iJ{pyP^^{-2igs)-^ -(2.3.3) 

The medium through which this highly off-shell parton is moving will now induce the 
additional radiation of a gluon that will carry away a fraction x of the parent parton’s 
momentum. The medium is modelled by GLV [11] as an ensemble of static scattering 
centres (that one can perhaps think of as a static, heavy parton), localized at say Xj = 
(zj, bj), that are all distributed with the same density 

p(q) = ^p{2). (2.3.4) 

The exact form of p{z) will have to be considered carefully, but for the time being it is 
only necessary that it be normalized such that / dzp{z) = 1. Each scattering centre is 
modelled as a Debye screened (or color screened Yukawa) potential with Fourier and color 
structure given by 


= 2n5{q^)v{q^,q^Je-^^-^ 


^Ta^{R)0Ta„{n), 


(2.3.5) 


where 




v{qn,Qn) = 

d-TTOs 

iW+Wn 


dvra^ 


4:7ras 

{Qn - ipn){q^ + iPn)' 


(2.3.6) 
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and ^,n is defined as 


fI = fJ-lx = + qn- 


(2.3.7) 


The effective small transverse momentum, differential elastic cross section is then given 
by 




CRC2{n) 4a2 


dA + 
CRC2{n) |u(q)p 
dA (27r)2 


(2.3.8) 


We are doing an opacity expansion (an expansion in the number of scatterings that occur) 
and therefore, to first order in opacity (scattering only once). 


d^N 


1 

— Tr \,y^Q T .y^\ T ^y ^2 T • • 
dt 


2 d^p 
(27r)32|p 


d^No + TV [{JKiJKl} + 23f?Tr{.^2^o } 

CLT 


d^p 

(27r)32|p| 


+ . . . , 


(2.3.9) 


where is the sum of all diagrams with one interaction with a scattering centre and 
^2 the sum of all diagrams with two interactions with a (single) scattering centre. First, 
to obtain the unperturbed inclusive distribution of jets in the wave packet, take the color 
trace of \^o\‘^ multiplied by the invariant one particle phase space: 

d^iVj = TV |.^olV 

(27 

= \J{p)\"dR- 

(2 



the spectrum of which can now be extracted as in [10] 


2.F(27r)3 2a;(27r)3 


d^Njd^Nj,^\ 


(2.3.11) 


Using equations (2.3.2), (2.3.10) and (2.3.11), one can obtain the radiation spectrum for 
the vacuum radiation 


diV^_ divj°^ CrOs 

dxd^k vr^ (k^ + 


(2.3.12) 
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In a similar fashion one can obtain the radiation spectrum for the first order in opacity 
energy loss due to final state medium interactions: 


1 


( 2 . 3 . 13 ) 


for which the energy loss is given by 


d^i 


d^k 






(27r)32’ 


( 2 . 3 . 14 ) 



Chapter 3 


Derivation of Energy Loss Formula 


In this chapter we will derive an expression for the first order in opacity energy loss of 
a high momentum particle as it propagates through a plasma of thickness L, mean free 
path Xmfp = l/(<7/o) and with color electric fields screened on a scale /i. In order to do so, 
one must calculate the invariant matrix elements of each relevant diagram and perform 
a number of manipulations in order to compute equation (2.3.14). Here we present the 
computed invariant matrix elements of the relevant diagrams and then sketch an outline 
of the process whereby equation 2.3.14 is computed. The details of the calculation are 
presented in the Appendices, sections B and C. 

We will first compute the invariant matrix elements of the Feynman diagrams associated 
with the leading order terms of the Dyson series expansion. The amplitudes are calculated 
by computing the / dqf integrals, closing the contour below the real axis (because zi > zq, 
meaning that zi—zq > 0 and so the exponents that appear will only converge if the contour 
is closed below the real line). To perform these integrals, one must find the poles, their 
residues and therefore the integral. I present here only the results of the computation of 
the amplitudes, please see Appendix B for details of the results below. 


Ml AO .ffk.c 
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'I ^> 
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k,c Mi_o.i 


S,Zo Z Zi 
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^ P 

‘^(r(r(rGf!r(r(r(r(r' , 

I - k,( 


Figure 3.0.1: Three “direct” terms -^1,0.1 contribute to the soft gluon 

radiation amplitude to first order in opacity L/A oc aei/A± 
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The matrix elements for the diagrams in Figure 3.0.1 are calculated in section B.1.1, B.1.2 
and B.1.2 to be 


^1,0,0 

e • k 




. (2vr)2 


d^Cll „-qgbi 


v{0,qi)x 


.^ 1 , 1,0 =J{p)EP^°{-i) f ^v{0,ci,)e-^^^-^^{-2igs)x 


(2vr) 


k • e 




^i{u}o+‘^m)izi-zo) _ f^-VlDi-z^o) 


Ctti 


^ 1 , 0,1 -J{p)E^P+’^>^{-i) J j^v{0,q,)e-^^^-^^2igsX 


e- (k-qi) 


(k — q;^)^ + + rrig 


X e 


i(a;o+(I>m)(^l-2o) _ D(‘^0-‘^i){zi-zo) 


c,ai]Tc 


ai 


(3.0.1) 


(3.0.2) 


(3.0.3) 


Notice that it is only the diagram that contributes a short separation distance 

correction. 

In order to properly perform the Dyson series expansion, we must consider diagrams with 
two interactions with one scattering centre, described in DGLV as two scattering centres in 
the ‘contact limit’ for the first order in opacity calculation. There are 2”’ — 1 diagrams for 
n scatterings, but we will find that not all 7 contribute. The diagrams that do contribute 
are shown in Figure 3.0.2, with the contact limit taken in cases where both the contact 
limit and well separated cases are shown. 

The amplitudes for the contact limits of the relevant diagrams are calculated in Sections 
B.2, B.3.1, B.3.2 and B.4.1 and given by equations (B.2.30), (B.3.29), (B.3.62), (B.4.22) 
and (B.4.23). Surprisingly, only the ^ 2 , 2,0 diagram contributes to the short separation 
distance generalization. 
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M 


2,0,3 

Z Zj 


^2,0,3 


Zt 


Zq Z Zj 


Zl 




^ p 

4 4 k,c 




P 

‘^(r6'f^6'()f6'6'<T(rf^^6'(r6'' 

k,c 


qi,ai 


(^2,3? 


qi,ai 



q2,a2 


(a) ^ 2 , 0,3 “direct” contributes to second order in opacity while 3 contributes to first order 

in opacity. 




k,c 


M 2 ,0,0 

Zo Z2 

^-1- 


z I 


% 




p 


M 2 ,0,0 

Zo ^Zi 

^ • i 


k,c 


- ^ § - 

qi,ai^ q2^a2 



qi,ai 


qT,a2 


(b) ^ 2 , 0,0 graphs for the well separated case as well as the contact limit. 


k,c 
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2,2,0 

Zo 


qi,ai 


Zl 
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# 




P 


§) - 

q2,a2 



(c) ^ 2 , 2.0 diagrams in both the well separated and contact limit case. 
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2,0,1 


®^Zo Z 




Zl 
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i p 
k,c 


qbai 



q2,a2 
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2,0,2 

^Zo z 






qbai 



k,c 


q2,a2 


(d) Diagrams .^ 2 , 0,1 and ^ 2 , 0,2 showing only the special contact case. 


Figure 3.0.2; Two scattering centre diagrams, taken from [10] 
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^2,0,3 


j{p+k)xo(_^^2 [ f 




(27r) 


1 


. (27r)2' 

e • (k - qi - q 2 ) 


X [[c,a2],ai](Ta2raO{e*(‘^°+“'")("i-"o) 

f d^^2 

(27r)2i (27r)2 


=. J(p)e«»+‘>™ / / |^•'(0. <ll)”(0. q2)e-‘'‘“+‘‘'>-‘’‘: 




02^01 


‘^2^2,0 ^J{p)s 


i(p+k)xo f d qi I" d q 2 -i(q,q,).b] 

J ( 27 r )2 J ( 27 r )2 


X ca2ai(r„,r,J 2 m2 + k2 + 


x(0,qi)'(;(0,qi)x 


/ 11 ^ p~U2(^l~^o) \' 

J(uJo+(llm)(zi-Zo) _|_ g-Ail(^l-Zo) ( _ \ 

V 2(/ii + P 2 ) ) 


^ 2 ^, 0,1 


=J(p)e*(P+^)^°(-i)2 
X e-^^^'^^2ig, 


d^^.1 „-iq,.bi ./n „ N f d‘^(l2 


(27r)2 

e- (k-qi) 


e-*‘b-^ix(0,qi)x(0,q2) 


(2vr) 


(k-qi)2 + + m^ 


gi(cJO+<im)(^l-2o) _ ^i{uJo-‘^l)izi-Zo) 


02 [c, ai](raira 2 ) 


^2^,0,2 


f ^^e-*‘b'^ix(0,qi)x(0,q2) f £^x 


X e 


-iq2-b2 


2i5s 


(2vr)2 
e- (k-qi) 


(2vr)2 


(k — q^)2 + M^x^ + m^ 

ai[c, 02 ] (Taj Taj. 


gi(tJO+“m)(2l-2o) _ gj(wo-6i;i)(2:i-2:o) 


(3.0.4) 


(3.0.5) 


(3.0.6) 


(3.0.7) 


(3.0.8) 


Here we will calculate the first order radiative energy loss. According to equation (10) 
[10], the starting point is the following quantity 

= (i T.-{|.//.P) + (3.0.9) 

is the sum of all diagrams with one scattering centre and ^2 the sum of all diagrams 
with two scattering centres. Therefore, in order to calculate this quantity, we must sum 
the relevant diagrams first. Once the summation is done, we square the amplitude (and 
take the real part of the trace), average over initial and sum over final states. This means 
that there is an averaging over impact parameter (introducing a factor of 1/Ax) and a 
sum over scattering centres. 
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In this section I outline the process whereby one obtains an expression for the first order 
(in opacity) energy loss by using equation (2.3.14). 

To simplify the process below, consider the following short hand: 

e • k 

_ £• (k-qi) 

(k — q]^)2 + 

= (wo - ojm){zi - zo) 

= {uo - uji){zi - zo) (3.0.10) 

= Woi - WOm = (wo - UJl){zi - Zo) - (wq - Wm)(^l “ ^^o) 

= -(wi- LVm){zi - Zo) 

= 

2 

= Tr(c^o^ — caca). 

Using the results from sections B.1.1, B.1.3 and B.1.2, given in equations (B.1.8), (B.1.25) 
and (B.1.15) and performing a rearrangement of terms so as to group like phases together, 
we obtain the sum of the one scattering centre diagrams. 


fk 

fq 

WOm 

WOl 

Wlm 

as 

a 


dx 


= N\J{p)\\Agi) 


2,^ 2x 1 f d\l C2{T) 


At J ( 27 r )2 


X < 4a/^(l - cos wom) + 2a/|(l - cos Wom) 


- 4a/fc/q(l - coswim) + a/,j/fc2 COS wqi + Trc^a^/|+ 


+e 


—^Az 


fl Tr (cos uom - 1) - Tr c^a^/| cos ujom 


+ fkfqa{cOS LOOm “ COS Wqi) 


+ ^/|TVc2a2e-2/^^^|. 


(3.0.11) 


Consider now the results in equations (B.2.30), (B.3.29), (B.3.62), (B.4.22) and (B.4.23). 
They can be summed, squared and averaged in a similar fashion to obtain 
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N\J{pf{4gl)L I 


d^<i, C2{T) 

(27r)2 dA 


X 


fl{2acosojQrn — 2a — Tr c^a^) 


+ 2a/fc/g(cOS Wom - coswoi) 
+ Tr cos wom - 


(3.0.12) 


It now remains simply to add the contributions from equations (C.1.12) and (C.2.3), 
which, after a number of cancellations and substitution into equation (3.0.9), finally gives 


. ^(1) ^ CnasLE f 
^ind 






ttA 


,g J TT (/i2 + q2)2 47r 


dzip{zi )X 


- ‘^fqifk - fg){l - COSUJim) 
p-uDi-zo) 


+ 


“ If) 


+ fkfgi COS OjQm “ COS Wqi) 


(3.0.13) 


where p{zi) is the distribution of the scattering centres which we will choose at a later 
stage (see equation (4.0.2)) and \g is the gluon mean free path. The second and third 
lines of equation (3.0.13) are then the short separation distance correction terms to the 
DGLV result. A cursory glance will reveal that our result behaves as can be reasonably 
expected: the dimensions of the correction terms are correct and the correction terms are 
suppressed exponentially in the limit that either Az —)■ oo or —>■ oo. That is, equation 
(3.0.13) reduces to the DGLV result in the large separation distance limit. Furthermore, 
as expected, equation (3.0.13) also goes to zero as Az —)• 0 due to destructive LPM 
interference. 

However, a surprising and wholly unpredictable cancellation occurs here; coincidentally, 
for the very same reason that many short separation distance correction terms were sup¬ 
pressed at the amplitude level, we see again that the correction terms are suppressed under 
the assumption that the formation length be longer than the Debye screening length of 
the scattering centres. It is noteworthy that, if one continues to assume that /r S> Wj, a 
term arising from the sum over single interaction diagrams cancels exactly with a term 
from the double interaction diagrams, leading to a correction term equal to 0 under the 
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assumption that /i » cjj. To relax the assumption that fj, ^ uji and rederive all expres¬ 
sions, will constitute a monumental task as the assumption led to simplifications already 
in the DGLV calculations and caused a massive reduction in the number of terms that 
need to be considered at the amplitude level in the present calculation. 

Lastly, note the breakdown of color triviality due to the CrICa term in the second line of 
equation (3.0.13), suggesting that, should fi ^ uji not hold (causing the retention of the 
correction terms), one will see a disintegration of the purely gluon final state interpretation 
in [12]. 

For completeness, the full energy loss formula (that remains to be integrated), is 


(1) _ CnasLE f d^gi 






7rX„ 


TT [gfl + q^)2 dvr 


dAzp{Az) X 


2(1 - cos {(wi -I- ojrn)^z}) 

(k — -|- 771,2 

(k-qi)-k (k-qi)2 


m2 -I- k^ -I- x 2M2 (k — q]^)2 -|- + m? 


+ le--‘A«){( 


m2 -|- k^ -|- x2M2 


X ( 1 - “ cos{(a;o - C 0 m)^z} 




k- (k-qi) 


(k^ -|- m2 -|- x^M^) ((k — q]^)2 -|- M‘^x‘^ + m'i) 


X (cos{(aio ~ L}m)Az} — cos{(a;o — wi)Az}) 


(3.0.14) 













Chapter 4 


Results and Conclusions 


In this section I present the results from a numerical analysis of equation (3.0.13), an 
improved energy loss formula, building on the DGLV result to include the effects of short 
separation distances between production and scattering. The computation in the previous 
sections of the energy loss was motivated by a desire to better understand energy loss in 
small systems, prompting us to consider possible modifications to existing energy loss 
formulae that could take small system sizes into account. Our naive approach was to 
equate the concept of a small system (that is, small L) to the idea of small separation 
distances (small Az), justified by the intuition that if energy loss occurs in a short, 
thin medium, the distance between production and scattering of the hard parton must 
necessarily be small as well. That is to say that, in order to generalize DGLV, we changed 
the length scale of the problem, taking: 

- < Am/p <C L ^ - <C L (4.0.1) 


where \mfp is the mean free path in the QGP. Therefore, while the calculation was 
performed with a focus on separation distances, the separation distance is eventually 
integrated over and one must look at the effect that the correction terms have on the 
length dependence of the energy loss formula. 

Presented here are plots of the numerical evaluation of equation (3.0.13) after analytically 
performing the integral over scattering centres. In computing the dzi integral, we assumed 
a distribution of scattering centres p{zi) such that 

foo POO f Az) 2 

J dAzp{Az) = J dAzexp I — 2-^| —. (4.0.2) 
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E (GeV) 

Figure 4.1.1; Energy loss per unit energy as a function of the parton’s initial energy 
for Leading Order (dashed) and Leading Order with NLO corrections (solid) for bottom 

and charm quarks. 

Analyses of light quark quenching at 130 GeV AuAu [35] suggest that the effective static 
opacity of the plasma can reasonably be fixed at L/\mfp ~ 4. Note also that the finite 
masses of the quarks considered in [10] shield the singularity at k — )• 0 which allows for 
numerical integration with momentum cut-off at 0. The numerics employ the following 
values: /r = 0.5 GeV, L = 4 fm where a constant for the length of the system was 
required, £' = 10 GeV oi E = 100 GeV where a constant parton energy was required, 
^mfp = 1 fill) Cr = 4/3, Ca = 3, = 0.3, mass of the charm quark was assumed 

to be 1.3 GeV while that of the bottom qark 4.75 GeV and lastly the QGD analogue 
of the Ter-Mikayelian plasmon effect was taken into account by setting the mass of the 
gluon as /r/-v/2 [36, 37]. All values were chosen to be consistent with the DGLV plots 
[10]. As in [10], kinematic upper limits were used for the momentum integrals such that 
0 < /c < 2x(l — x)E, ensuring collinearity, and 0 < g < ^/3E^, due to finite kinematics. 
The fraction of momentum carried away by the radiated gluon, x is integrated over from 
0 to 1. 

4.1 Numerical Results 

Figure 4.1.1 shows results from numerically integrating equation (3.0.14) as a function of 
the initial energy of the hard parton. We compare the leading order energy loss of charm 
(because it is an experimentally clear probe [38]) and bottom quarks to their energy 
loss when small separation distances are considered. In the figure captions, NLO refers 
to Next to Leading Order in separation distance (as compared to the Debye screening 
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Figure 4.1.2; Energy loss per unit energy as a function of the static thickness of the 
medium for Leading Order (dashed) and Leading Order with NLO corrections (solid) for 

bottom and charm quarks at 10 GeV. 

length). Note that the results show a decrease in energy loss due to a negative correction 
to the LO result. 

The plot in Figure 4.1.1 differs slightly from its counterpart in [10] due to somewhat 
different kinematic limits in the numerical integration. 

Figure 4.1.2 shows the differential energy loss of charm and bottom quarks as a function 
of the static thickness of the medium for both the LO and the NLO (LO with small 
separation distance correction). Again we see here a reduction in energy loss to the point 
that we expect to see a small amount of particle enhancement for systems less than a 
femtometer across. A striking feature of the thickness (L) dependence of the correction 
term is that the effect of including short separation distances does not seem to decrease 
with increasing system size, suggesting that the effect should be taken into account even 
when performing calculations in the large system size approximation. The fact that the 
NLO correction applies even to large systems is not a complete surprise as the dzi integral 
is performed over all zi, from 0 to oo. 

However, when one considers the same plot for 100 GeV partons, shown in Figure 4.1.3, 
the negativity of the energy loss becomes so large that one must revisit the premises of the 
calculation. The numerics seem to suggest no suppression whatsoever, even for sizeable 
systems - a postdiction that simply does not describe the plethora of available data for 
particle suppression in AA [4, 6, 39]. 
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Figure 4.1.3; Energy Loss per unit energy as a function of the static thickness of the 
medium for Leading Order (dashed) and Leading Order with NLO corrections (solid) for 
bottom and charm quarks at 100 GeV. 

4.2 Discussion and Conclusion 

In order to interpret these results, it is important to highlight a key feature of the calcula¬ 
tion of the correction terms, which involves the assumption that the formation time of the 
radiated gluon is much larger than the Debye screened length of the scattering centres. 

In the large separation distance limit that was employed by DGLV [10], the effects of the 
large formation time assumption were fairly subdued due to the much larger and more 
immediately evident exp{—Az/r} exponential suppression. However, when the calculation 
is performed under the assumption that the separation distances are on the order of the 
Debye screened length of the scattering centres, one must deal explicitly with a regime 
of /S.Z values that are small enough that the formation length becomes the largest length 
scale^ of the problem. The small separation distance regime necessarily enhances the 
effect of the large formation length assumption. 

The large formation length assumption introduced a considerable calculational simplifica¬ 
tion, leading to a surprising and non-trivial cancellation effect. The suppression of a large 
number of NLO terms appeared already at the amplitude level of the calculation, leav¬ 
ing only two amplitudes, one single scattering and one double scattering, with non-zero 
NLO terms - even under the assumption that fj,i 3> wi. However, when performing the 
ensemble average within the Dyson series, the //i S> cui assumption leads to a non-trivial 

^To be precise, when considering small separation distances, the mean free path is the largest length 
scale of the problem, but it is not relevant to the calculation because the mean free path is a fixed, 
calculable quantity that is considered to be much larger than the Debye length of the scattering centres, 
and therefore larger than the system size, which is the scale we work on. 
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cancellation, resulting finally in a null correction term to the end result. Numerically, if 
one does not assume that the large formation time approximation holds (as we have done 
here), the correction term contributes more than 100% as seen in Figure 4.1.3. Relaxing 
the /ii 3> oJi assumption would entail an immense complication of the calculation (even 
for large path lengths, the relaxation of the large formation time limit would include extra 
terms), but it is clear that the effect must be taken into account since, as will be elabo¬ 
rated on in the next paragraph, the small separation distance limit also affects the large 
system size calculation. 

Figure 4.1.2 clearly shows that the correction to the energy loss due to the inclusion of 
small separation distances remains sizable even for large systems. Upon closer inspection 
of the details of the calculation, the reader will notice that the integral over separation 
distances is in fact performed starting at Az = 0, meaning that small separation distances 
are included even if they are neglected in the calculation at the amplitude level. The fact 
that the separation distances are integrated over the entire range of distances suggests that 
the corrections are, in principle, important for large system sizes as well. The numerics 
presented here explicitly show the importance of the small separation distance correction 
terms, particularly at very high energies. 

The heavy ion physics community has stumbled across an excellent opportunity to scru¬ 
tinize existing interpretations of data due to surprising recent discoveries in what should 
have been calibrative exercises. The use of pA results as a probe of nature of the QCD 
coupling and the validity of competing theories relies heavily on a robust understanding 
of known processes. Such an understanding is contingent upon extensive analyses of the 
two most common signatures of the QGP, collective behaviour and particle suppression. 

We have endeavoured to catalyze the analysis of energy loss in pA by calculating the small 
separation distance corrections to the large path length static energy loss formulae and 
have found that: 

1. Under the assumption that the Debye screened length of the scattering centres is 
much larger than the collinearity of the radiated gluon, a remarkable cancellation 
occurs and there is no correction. In future calculations this assumption will have 
to be relaxed since it is not clear that pi S> ui holds throughout the momentum 
and collinearity ranges that are considered. 

2. If one relaxes the large formation length assumption after the amplitude stage of the 
calculation, one finds that, since the separation distance between production and 
first scattering is integrated over the entire separation range, the effect of including 
short separation distances is > 100% even for large system sizes and must be taken 
into account. 
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Useful Results 


A.l 

From [10], the amplitude for hard jet radiation to emit a transverse plasmon with mo¬ 
mentum, polarization and color (/c, e, c) is 


-h k)fj,e^{k)i Am {p + k)c 

In (A.1.1) we have used that xo(0,zo,0) is the jet production point inside the plasma. 
We also have that J varies slowly with p which allows for the above approximation since 
k ^p. 


A.2 Writing down the Matrix element 

In this section I look briefly at the ‘Feynman Rules’ for the invariant matrix elements we 
consider. Take for example the matrix element for .^ 1 , 0 , 0 ; which is given by 


^1,0.0 = / + + 

xi Am {p- qi + k)i Am {p - gi)(2p - qifV 

~ j(p + (A.2.1) 
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where 


h{p,k,qi,zi - zo) 


dqf ea(2p - 2g +/c)" 

27r {p — qi + ky — + ze ^ 


X 


{p — diY — 




(A.2.2) 


To see this, consider the first diagram in Figure (B.1.1). Equation (A.2.1) can be ‘read’ 
off by going from left to right and considering the following set of ‘Feynman rules’: 

• iJ{p + k — qi)e^^‘P^^~‘^k)^°^ which is the highly localized wave packet produced in free 
space, propagating with momentum {p + k — qi). 

• {igs)ea{2p — 2q + kyaicTa^, the radiation vertex with color factor . The form comes 
from ‘Peskin’s trick’ described in section A.9. 


• Am {p — (li + k) and Am ip — <?i)- These are the propagators of particles with mass 
M and with momentum {p — qi + k) and {p — qi) respectively. They have the form 


Am {p-qi + k) 


1 

{p — qi + ky — + ie 


(A.2.3) 


• V. The potential, the random color screened potential that is given in [10] 
as 

Vn = V{qn)U^"^" = 27r5((?0)u(q„)e-'^"-""ra„(i?) ® T^Jn), (A.2.4) 

where x„ is the location of the scattering centre and u(q„) = 47rQ!s/(q^ + 


A.3 Components of k, p and q 

Here I present the derivation of the components of the momentum (four) vectors fc, p 
and q. We are considering the diagram shown in Figure A.3.1. Recall that the invariant 
masses are such that p-p = and k-k = m^. We also know from conservation of energy 
that p + k = P. Suppose now that the gluon carries away a fraction of the initial quark 
momentum x. We then have, in light cone coordinates: 


k = [xP~^, k , kj; 
p= [y - x)P^,p~,-ky, 


(A.3.1) 

(A.3.2) 
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Figure A.3.1; A gluon with momentum k radiated from a quark with initial momentum 

P+. 


It remains only to determine k and p . Therefore 

p ■ p = p'^p~ — p • p 

= (1 — x)P^p~ — (k) — k 
= (1 - x)P+p- - (k)2 = 
m 2 + k2 
“ {l-x)P+ 


(A.3.3) 


k-k = k+k- - k-k 


= xP~^k —\sf = mfg 
m2 + k2 
-^oP^ 


(A.3.4) 


Since e • A: = 0 by the Ward identity, we can find the polarization vector for the emitted 
gluon by choosing = 0; 


e • k 


- (e+/c — e k~^) — e • k = 0 

_ 2e • k 

^ ^ k+ 

2e ■ k 

xP+ 


(A.3.5) 


We therefore have that 


k 


xP+, 


nig + k^ 

xP+ 


k 


p = 

e = 


{l-x)P+, 


0 , 


2ek 

xP+ 


e 


M2 + k2 
{l-x)P+ 



(A.3.6) 

(A.3.7) 

(A.3.8) 
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A.4 Recurring dot Products 

In this section I present results for three dot products that appear regularly. These are 
the dot products of four veetors. 


A.4.1 2p ■ k 


Using the definition of the dot product in light cone coordinates given in Section 2.1 


2p-k = p~^k +p /c'*'— 2p • k 

^(1 — x)P~^{ml + k^) 


xP+ 


9 

m? + k^ 

-+ (M2 + k2)x-2k2, {x<.1^{1-x)Ak A) 

X 

1 


+ 


{l-x)P+ _ 


(M^ + k2)(xP+) 


-2(-k)(k) 


= - + k^ + x^M^ + k^x^ - 2k^xP 

X ^ 


[m^g + k^ + x^M^], 


(k2(l + X + x2) «k2) 


(A.4.1) 


A.4.2 2p ■ q 


Dot products with q ned to be started off slightly differently due to a slight simplification. 
We use the definitions set out in Section 2.1 to obtain. 


2p-g = 2(j>V-pV-p-q), (g° = o) 


= 2 - -P ) 


- 2p • q 


(1 - x)P+ - 


2 n 


m 2 +k 

(1 -x)P+ 


g" - 2p • q 


-P+g"" + 2k • q 


(A.4.2) 


In the last line we have used that 


^ - P ) 


IL M2 + k2l 

2 L (1 

ir + M2 + k2i 

- x)P+_ 

^P+ 

2 

2 ^ P+ 


(A.4.3) 
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A.4.3 2k-q 

In a similar fashion to the process that yields equation (A. 4. 2), we obtain 


2k ■ q = 2{^ - - k ■ q), 


= -2 


xp+ _ ml + 

xP+ 


(gO = 0) 
^ -2kq 


-xP+q^ - 2k 


(A.4.4) 


A.4.4 2p ■ e 

Recalling that e"*" = 0, we have the following, again relying on the definitions in Section 

2 . 1 . 


2pe = 2 


-(p+e“ - p • e 




(1 -x) 


X 

2e ■ k 

X 

2e ■ k 


2e • k + 2e • k 


— 2e • k + 2e • k 


(P = -k) 


(A.4.5) 


A.5 Poles from propagators 

A.5.1 AM(p-gi) 

Since it is the denominator of the propagator which has the form given in (A.2.3), we 
need to solve the following equation for qf (I drop the subscript 1 for simplicity since it 
will not cause any confusion): 


{P -q)^ - + ie = 0 


(A.5.1) 
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This implies that 

0 = — 2p ■ q + + ie 

= M^- 2{-p^q^ - P • q) + (-(g")(g") - q^) - 
= 2p^q^ + 2p • q - (q^)^ - q^ + ie 

^ P+q^+2p-q-{q^f-q^+ie (A.5.2) 


One can now do one of two things: One can solve the above for q^ directly, or one can use 
the quadratic formula. The latter option is not only more rigorous, but also allows for a 
simpler handling of cancellations that occur when the residues are calculated. Therefore, 
I will follow this method: 


Q = 


-P+T V(P+)^-4(-l)[.] 
2(-l) 

P+ 7(P+)2 + 4[.] 




= ^ ^ 1 , ^ 
2 ^ 2 V (P+)2 


P+ 

P+ 

~Y 

p+ 

~Y 


1T^1 + 

It (i + 
It fiT 


4[. 


( p +)2 

(p+)2^ 

2 


= 2p • q - q"* T ie 


(By the Binomial Theorem) 


[2p • q - q^ T i 


ie\ 


(A.5.3) 


Which gives us (as expected) two solutions. We ignore the solution from the plus sign 
because it gives Tie which lies in the upper half of the complex plane and is therefore not 
included as we close the contour below the real axis. The solution from the minus sign is 
then 


q^ = 


P+ 

~Y 

p+ 

~Y 

-1 


1 -1 - 


(p+)^ 


:(2p • q - q^ Tie) 


(P+)^ 


:(2p • q - q2 Ti 


^e 


= ^[^P-q-q Tie 

2k • q q2 

T — ie, 


P+ P+ 

(k T q)^ — k62 


p = —k , see equation A.3.7 


(A.5.4) 


— ze 


P+ 


(A.5.5) 
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Although this result differs slightly in form from the corresponding result in [10], the 
difference does not affect later results. 
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A.5.2 Am {p - Qi + k) 


As in the previous section, we need to solve the following equation for q\: 


0 = (p — + k)'^ — + ie 

= + (f‘ + 2p ■ k — 2p ■ q— 2k ■ q — + ie 

= + ml + q^ ^ --h P^q 


— + scP^^ + 21^-^ — + i 


le 


y , rN9 9 rUg + , 

= ^ - {q )^ - H- - -h P^q + ie 


(A.5.6) 


In the third and fourth lines we have used the fact that x <C 1 ^ xA <C A to neglect 
certain terms. We can now solve for q^ in much the same way as in the previous section: 


<1 = 


-P +±'/(P+P-(4)(-l)H 
2{-l) 


rUg + k^ + M^x^ g 
* = — -+ ie 

X 


P+ P+ I 4 r, 

r / 2 

+(7^1* 


(A.5.7) 


Which again has both a positive and a negative root and we will again ignore the root in 
the upper half of the complex plane. The root in the lower complex half-plane is then 


P+ 


/-X- 


(p+)2 


+ k^ + M^x^ 


X 


- q + ie 


r2^2 


1 /rrig + k +M^x „ 
^-" 


X 

2^2 


+ k^ + M'^x 
~ xP+ 

= £X^0 


- ^e, 


(q^ockX k2(l + x)f«k2) 


(A.5.8) 
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A.5.3 {k - qi) 

As in previous sections, we must solve the following equation for qf and we will drop the 
subscript 1. 

0 = {k — q)"^ — + ie 

= k^ — 2k ■ q + q^ — m?g + ie 

= y^g - 2 [0^ - k^q"" - k • q] + qf - + ie, (g° = 0) 

=-2 - ^(/c+- A;“)g^ - k • q + - q • q + ie 

= q^xP^ + 2k • q — — q^ + ie (A.5.9) 

Again we can use the quadratic formula to obtain two solutions: 



[*] = 2k • q — q^ + ie 


(A.5.10) 


Again we must choose the root that lies below the real axis: 



= ;f^(-2k-q + q) - ie 

= 7i^[(k-q)'-k'] 


= oji — uio — le 


(A.5.11) 
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A.6 Poles from the Potential 


The potential is given by equation (A.2.4) as 




The parts of this that we are concerned with are r’(q„), given by 




47ras 


47ras 

-( g ^) 2 - q 2+;.2 ^ 


Tn)) 


(A.6.2) 


for which we have to find the qf's that result in infinities. Using the definition /i? = = 

q? + we have that 


Tn + = q? + 

= (9D'-q? + /^? = o 

^ {qtf = q? - /r? 

2 2 2 2 
— qi “ qi “ /^i — ~ui 

qf = ii/ii- (A.6.3) 

However, both of these poles were neglected in the large separation distance limit. The 
positive pole is neglected because it lies in the upper half of the complex plane and the 
contour is closed below the real axis. The negative pole is neglected in some cases because 
its contribution will be accompanied by an exponential term of the form exp —i^{zi — zq) 
which is then exponentially suppressed because we deal with the well separated case in 
which ^{zi — zo) = /rA ^ 1. In the contact case related to the diagrams with more than 
one scattering centre this is not the case. 

In this treatment, where the size of the system is of the order of the Debye screening length, 
we will consider the pole in the lower half of the complex plane since its contribution will 
no longer be exponentially suppressed. 





Useful Results 


42 


A.7 Factorizations within residues 

A.7.1 {p-qf-M^ 


{p - q? 


— 2p • q + 

= -2(A®' - p^q^ - p • q) + (A®' - q"q" 

= 2p^q^ + 2p • q - (g")^ - 


= 2 




g^-2k.q-(g^)2-q 


qq), 


= (1 - x)P+g" - 2k • q - {q^f - q^ 


^ P+q^ + k^ -{k + qf 


= P+x 


g^ 

X 


+ ujo + Cj 



(g° = 0) 


A.7.2 (p-g + A;)2-M2 

(p — g + A:)^ — — 2p ■ q + 2p ■ k — 2q ■ k — 

= - q^ + mA - .2^^ 

+ l/x(m^ + k^ + x^M^') + xP^q^ +21^^ 

P“''g^ + l/x(m^ + k^ + x^M^), (x <C 1, P"*" » 1) 

= P^g^ + P'''(wo + 

= P^(g^ + Wo + Wm) 

In the third line we have made a number of cancellations based on approximations of 
large P"*" and small x. It is important to note that they can only be made as comparisons. 
That is, terms can only be neglected if they are much smaller than a similar term that is 
different only by a factor of x or P'*'. 
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A.7.3 {k — qy — mfg 


{k - qy - ny = j/ -2k ■ q + 

= —(xP~^q^ — 2k • q) — {q^y — q^, Section A.4.3 

= -xP+ + 2k • q - - q2 

Ri -xP+q^ - (k - q)^ + k^ 

= —xP'^q^ — xP^{u}o — wi) = —xP{q^ — oji + ojq) (A.7.1) 

A.7.4 v(q„) 


v{^n) 


47ra<j 


qn + M 

AtTOs 


■= v{qf,qn) 


(g^)2 + q2 + /i2 


?;(0,q„ 


We will encounter two kinds of residues: Ones that contain w’s and ones that contain 
—ill’s. The above result will hold in both cases. For the former, it holds due to the small 
X approximation. In the latter case it is identically true. We will also see, when dealing 
with diagrams that have two scattering centres, situations where the poles containing 
—ill are not suppressed. We will then see two cases: one where the z-component of one 
q vector is used along with the perpendicular component of another and a second case 
where the same q vector is involved. In the first case we have 




Anas 

-{qiy + 


Anas 


(giy + ^ii + Ui - q?’ 

Anas 

(g| + i//i)(gf -i/ii)’ 


il^l = Ui± = q? + 


(A.7.2) 


while in the second case we have that 
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(A.7.3) 


{Q2 +*/^2)(g2 -*^ 2 )’ 


A.8 Color Matrix results 


Using the relations set out in section 2.1, I prove here that 2i/“^^c6a 


LHS = = C2{G)5‘^cd 

= C2(G)cc 

= -G2(G)G2(r)l 



= G2(G)G2(r)i = -LHS 

+ 2ir^^cba = 0 


It is also worth noting that 


[A,B]* = {AB-BA)* 

= (AB)* - (BA)* 

= B*A* - A*B* 

= BA — AB Hermitian 
= [B,A] = -[A,B] 


We will also need to know that = a^c^. It is trivially true by change of variables, but 
can also be shown from the definitions: 
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ccaa = c(ac + = caca + if'^°‘^cba 

= (ac + i/“^6)ca + if^°‘^cba 
= acca + if‘^°'^bca + i/“^c5a 
= ac(ac + i/“^6)z/“'’6ca + if^‘^\bc + i/^“(i)a 
= acac + i^'^’^acb + i^^-^ca + if^^%ca - 
= a{ac + irh)c + ir\{bc + + 2ir%ca - 

= + if^^abc + ir\bc - f^^f’^^ad + 2if^^bca - 

= + 2if^^abc + 2if<^%ca - f^^f^'^ad - f'^^f^'^da 

= c?c^ + - 2f^^f^^da = c?c^ 

Also note that, by change of variables, acac = caca. 

A.9 Peskin’s Trick 

We will be making use of ‘Peskin’s trick’ as it is called in [40] . It pertains to the Feynman 

rules of diagrams with bremstrahlung and appears in Peskin on page 183.: 

)7''e* (/ + m) = u{ph''e*^{Ypl + m) 

= u{p) [l^l'^elp^ + 7^e* m] 

= u{p)[{ — 7 ^ 7 ^ + 2g^'')e*^p^ + 7 ^e*m], (Antiommutation relations) 
= u{p)2e*^pl - Yl^u{p)elp^ + u{p)-i^elm 
= u{p') [2e*p' + {-f + m )7 • e*] 

= u{p)2e ■ p', (Dirac Equation u{p'){f> — m) = 0 ) (A.9.1) 

The conjugate relation follows exactly the same reasoning. 

A. 10 Major Simplification Effect 


The vast majority of contributions from previously exponentially suppressed terms are 
still suppressed but by a different factor. I discuss here that effect: 
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This entire calculation (even the original) worked under the assumption that uji <C /xi- 
That is to say that 

k2 

This in turn implies that 

11 ^ 
— > — A.10.1 

This fact is used very often in the calculation and is true for all uJi and /Uj. 



Appendix B 


Detailed 


calculation of Amplitudes 


In this chapter I present the details of the calculation of the amplitudes of the relevant 
Feynman Diagrams. These calculations closely follow the Appendices of [10], with the 
exception of additional terms that are considered in the small separation distance limit. 


B.l ^1,0,0? -^1,1,0? ^1,0,1 

The basis results appear in [10] in Appendix B. The diagrams are shown in Figure (B.1.1). 
In what follows, I will often drop the subscript 1 when using q {qi = q). 

The general idea behind what is done here is to calculate the f dqf integrals by closing 
the contour below the real axis (because zi > zq, meaning that zi — zq > 0 and so the 
exponents that appear will only converge if the contour is closed below the real line). To 
be able to do this, one must find the poles, their residues and hence the integral. Many 
of the results are repeated and are therefore presented in Appendix A. 


Ml.0,0 

/ 


I qi,ai 


M 


1.1.0 












k,c M 




Z 

I- 

I qi’^i 


1.0,1 


^Zq Z 


Zl 


^ P 

M - k,( 

I qh^i 


Figure B.1.1: Three “direct” terms contribute to the soft gluon 

radiation amplitude to first order in opacity L/A oc 0 ^. 1 / A±_[\Q]. 
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B.1.1 Computation of 

The matrix element for ^ 1 , 0,0 is given by section A.2 as 


^ 1 , 0 , 0 = / ^iJ{p + k-q,)e^(P+’^-‘^^>°{igs)ea{ 2 p- 2 q + krx 

X iAM{p -qi + k)iAM{p - qi){2p - qifV{ql)e"'^^^^Ta^aic 

^J{p + k)e^^P+^^^%-igsaicTa,)2E J (B.1.1) 


where 


^i(p, -zo) 


dql ea{2p - 2g + A:)" 

27r {p — qi + k-Y — + ie^ 


X 


(p “ 91 )^ “ 


^'(qi 




We can simplify Ii slightly through noting that 


ea{2p — 2q + k)°‘ « 2pe 


2e • k 

X 


(see Section A.4.5) 


(B.1.2) 


To solve Ii, note that there are three poles; one from each propagator and one from 
the potential. These poles are calculated in Sections A.5.1 and A.5.2. (The third pole, 
calculated in Section A.6 is no longer neglected for the reasons explained there). The 
poles are then, from the propagators 

OJq 

(k + q)2 - k2 
-- 

ipi 


z(l) 

<ll = 

z(2) 

= 


z(3) 

= 


We must now calculate the residues due to these poles, taking the limit as e —0. The 
residue for the pole at q^^^^ = —uq — Cjm — ie is given by: 


Res[g^^^^ = -Wo -iom- ie] 


= lim 

^ 2 


( 1 ) 


^{qf-ql 




v{e[i)e 


{p — q + kY — + ie{p — qiY — + ie 


(B.1.3) 


Each part of this can be evaluated individually, by sections A.7.2, A.7.1 and A.7.4 
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2e-k 

X P+q^ + k2 - (k + q)2 

x(P+)2(a;o + cDm) 


In a similar way we can calculate the residue for the pole at which is 

given by 


Res[/i,( 7 ^*'^^] = lim 


z 2 ( 1 ) 






2e • k 


{p — q + ky — + ie{p — gi)^ — + ie 

exp{ - f (l"+q|-'^ (^i 2 :o)}^;( 0 ,q) 


(k>q)^ 

P+ 


2e-k exp{z ^^+^^ ^ (zi - 2 :o)}u( 0 ,q) 
x(P+)2 (wo+Wm) 

2e • k u(0,q) 
x(P+)2 (wo +Wm) 


(B.1.5) 


We now consider a residue that was not calculated in [10] because it is exponentially 
suppressed in the large separation distance approximation. The residue due to the pole 
Sit qf- = —^/il, IS given by 


T 5 I z(3) • ri r c[lJrdqIf{47ras) 

Res[gi' = = hm , ^ ^ -r- 


2e • k 


ialA-rpffiqf-ifil) X 

^-iqf(zi-zo) 


(p — q + kY — +ie{p — qiY — + ie 

( 47 ras) 2e • k 1 g-Mi(zi- 2 :o) 

-2i/ii X P+(-i/ii+Wo + Wm) R+(-i/ii) Tk^ - (k + q)2 

(47ra,)2e-k 1 1 

—2ip,i X {P~^Y 
47ras 2e • 


2x(P+)2 i(^3) 

47rao 2e • k 


x(P+)2 /i^ 




(B.1.6) 


These results can be multiplied together to find Ii. Also note that the small x assumption 
reduces the exponent in equation (B.1.5) to 0 in comparison to the exponent in equation 
(B.1.4). Since the integral is from —00 to 00 and we close in the lower half of the complex 
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plane, the contour is negatively orientated. We therefore have 
—2Tri 


h = 


-27r 


— i 


Res 

2e • k ^(Ojq) 
X{P+Y {(jjQ+Cbm) 


( 1 - 


^i{ujQ+u} 

m ){zi-zo 






x(P+)2 flf 
2e • k dvros 

y(P +)2 /rf 


gi(a;o+i^m)( 2 i- 2 o) _ g-^ti( 2 : 1 - 20 ) 
- i - 


(cjQ + Wn 


IJ-l 


(B.1.7) 


where we have used the definitions of a;, E, and E~^ in section 2.1. Here we encounter for 
the first time an effect that will cancel many of the contributions in the small separation 
distance limit. The effect is discussed more fully in Section A.10. Since wq <C /j-i, the 
second term in the large brackets - which is also the ‘new’ term - is suppressed in relation 
to the first. This means that Ji remains unchanged under the small separation distance 
assumption and therefore so does the amplitude: We now have the matrix element 


^ 1 . 0 , 0 -J{p)E^P^’^^^^{-igs)aicTa, 2 X J 

e • k 


„-q,. 

(27r) 


e 


u(0,qi) 


^ (k^ + 777,2 + M2 x2) 

=J{p)EiP+>^')-o^2ig,)aicTa,{-i) J 


d'^^1 - 


e • k 


k^ + 7772 + M2 x2 


(2vr)2 

|‘gi(wo+Wm) _ 


qi-bi 


77(0, qi)x 


B.1.2 Computation of ^110 
The matrix element for is 

^ 1 , 1 , 0 = J ^iJ{p + k-qi)E^P+'^-^^>°{ 2 p + 2 k-q,fx 

X iAM{p + k - qi)iAM{p + k){igs)ea{2p + k)°‘V{qi)e'''^^^^Taj^cai 

^J{p + k)E<^P+^^^°{-ig,Ta,ca,){2E + 2u;) J x 


1 


1 


{p + k — g'i)2 — m 2 + ie (p + k^ — M2 + ie 
The same simplifications apply as before. We have that 

2e • k 


(27r)’ 

-ea{2p + kY 


(B.1.8) 


ea{2p + A:)" « 2pe 


X 


(see Section A.4. 5) 





















Detailed ealeulation of Amplitudes 


51 


while the second propagator does not depend on qi and simplifies as per equations 
(A.4.1)and (A.4.5) to give 


ea{2p + ky ^ 2p • e 

{p +ky - M‘^ + ie ^ + 2 p -k - 

2p ■ e 

k"^ + 2p ■ k 

2k•e 1 

X 

_ 2k - e 

772,2 _|_ Jj-2 _|_ 2-2^2 


(B.1.9) 


We also know that, since uj E, the simplification (2E + 2uj) ~ 2E applies. We therefore 
have that 


.^ 1 , 1,0 = J{p)EP^°{-igs)Ta,cai4E 


k • e 


m2 + k^ + 


w q 

rC 


(27r)2 


*‘^''’'/ 2 (p,A:,q, 2 ;i - zo), 
(B.1.10) 


where 

hip.k.a, -z„) = l + (B-i-ii) 

There are only two poles: One from the remaining propagator and one from the potential. 
Again we now retain here the pole from the potential that had previously been suppressed 
exponentially. The pole from the propagator is calculated in section A.5.2 to be ~ 
—Wo — Wm — Taking the limit as e —)■ 0 and using results from sections A.7.2 and A.7.4, 
the Residue due to the pole at is 


Res [l 2 ,q^^^^] = lim 


(p-q + kf-AP + ie'^^’ 


q^ + Up + Um i(ojn+CA^)(z, -Zn) 


P~^{q^ + Wo + OJrn) 


2 ;( 0 ,q), 


1 

p+ 


m )(zi-zo). 


e""”“~"' 2 ;( 0 ,q). 


(B.1.12) 


For the same reasons as for the previous residue, we have the residue due to the pole at 
= -igi 

Res\a^^^h= lim [ 

^ - ip.) P+{-qi + Wq + Wm) 

P+{-ipi){-2ipi) 


2P+(/xi)2 


(B.1.13) 






















Detailed caleulation of Amplitudes 


52 


We therefore have that I 2 is given by 


h{p,k,q,zi - zo) = 


P+ 


J{uJo+Cj 

m )(zi-2o) 




(4Tras)e ^ 1(^1 
2P+(/ii)2 


= . .^J47ra^ 

^ V?(2p^)+l 


gi(aJO+i^m)(2l-2o) _ _g-Ml( 21 - 2 o) 


(B.1.14) 


Where the negative sign is again due to the negative orientation of the contour (that 
closes in the lower half of the complex plane). The matrix element is therefore 


^ 1 , 1,0 ~ J{p)e^P^<^i-igs)Ta,cai 4 E 


k • e 




771,2 _|_ ]^2 _|_ j.2]^2 J (27r)2 


e *qrbi__x 


2p^ 


1 


i(aJO+“m)(2i-2o) _ _ -Ati(2i-zo) 

2 


77 ( 0 , qi) 




J (27r)2 

k • e 


77 ( 0 , qi)e 


-iqi-bii 


m2 + + x2M2 


{-2igs)x 


1 


gi(a;o+tjm)(2i-2o) _ _g-7il(21-2o) 


Ta^cai (B.1.15) 


Here we see then for the first time the retention of a term in the small separation distance 
treatment that is exponentially suppressed in the large separation distance approximation: 
Notice that the second term in square brackets is exponentially suppressed in the large 
separation distance limit, meaning that B.1.15 reduces to the large path length expression 
in the limit that l//ii <C (zi — zq). 


B.1.3 Computation of 

In a similar fashion, the matrix element for .^ 1 , 0,1 is given by 

iAM{p + k - qi){-i)Amg{k - qi) 

^J{p + k)P(P+>^^^o[^,ai]TaA-i) J 

X 2^56 • (k - qi)2£'l3, (B.1.16) 

where 

h{p,k,qi,zi - Zo) = J ^v{ql,q^)AM{p + k - qi)Amgik - (B.1.17) 


To solve Is, we again consider the poles from the potential and the two propagators. We 
now consider all three poles, including the one from the potential. These are, as per 
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sections A.5.2 and A.5.3, at (dropping the subscript) 


qAD 

qA2) 

qA3) 


— UJq i^ 

= —Wo + wi — ie 
= -itH 


(B.1.18) 


We need to calculate the residues due to each pole. Taking e —0, the residue due to the 
pole at 


Res = lim 

gZ^gZ(l) 


{q^ — — Cjm) e ^o) 


{p + k — qY — + ie {k — qY — w-g + i^ 

Again we can simplify using sections A.7.3, A.7.4 and A.7.2 so we have that 


(B.1.19) 


Res = -Wo - w^] 


xP+_(gl^h-woH=T^^ xP+(-wo - Wm + Wo - Wi) 


_ 1 '^(0;^) -i(-u!fy-a,^)(zi-zo) 

X{P+Y + Wl 


(B.1.20) 


We can use exactly the same methodology to find the residue due to the pole at = 
—Wo + Wl — ie. 


Res 


lim 

qz^qz{2) 


(g2 — g,z(2))^(— ^0 + Wl) e Zq) 

{p + k — qY — + ie {k — qY — Tn'^ + ie_ 


(B.1.21) 


With the same simplifications that were used for the previous pole, we have 


Res = —Wo + wi] 


P+(-wo + wi + wq + Wm) 

1 'U{^0 -^^o) 

x(P+)2 cD 

m + Wl 


(B.1.22) 
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The third residue is given by 


= lim 


z V 2(3) 


iaf-A-ifa) {A-iras 


1 


{qf - ip.) {p + k - qi)'^ - + ie 

^-iq)l{zi-zo) 


{k — qi)'^ — + ie 


1 


^-fil{zi-zo) 


(47ras) ____ 

{-2ipi) P+{-ipi + Wo + Wm) {P~^x){-ipi - wi + wq) 

dvra, 1 

{-2ipi) (x(P+)2) (-i/ri)2 


1 dvroc 


2(-i/ii)3x(P+)2 


„-/^l(zi-^o) 


(B.1.23) 


The integral I 3 is then (including a negative sign for the negative orientation of the 
contour) 


h{p,k,qi,zi “ ^ 0 ) ~ 


dvra^ m )(zi-zo) _ gi(cJo-tJi)(zi-zo)^ 

x(P+)2/xf(a;i + Wm) 


1 dvra^ e ^ 1(^1 ^ 0 ) 1 

+ 2^,;f (i/ri) x(P+)2 


(-i)(47rQ;s) 

^x(P+)2(^2) 


oi([»;o+<im)(zi-2o) _ g*(aJo-t*Ji)(2l-2o) g-/ii(zi-zo) 

-+ 


(wi + U)m) 


{ 2 ipi) 


(B.1.24) 


We see again here the effect described in A. 10, resulting in the second term in the large 
brackets being suppressed in relation to the first, meaning that this integral remains 
unchanged under the small system assumption. The entire amplitude will also remain 
unchanged, we therefore have the Matrix element 




d 2 (ii 


u(0,qi)e 


e- (k-qi) 


(k — q;^)2 + + m2 


X e 


(< 


,i(a;o+tjm)(2l-2o) _ g*(<.^0-‘<.’i)(2i-2o) 


[c, aijTc 


ai 


(B.1.25) 
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M 


2,0,3 

Zq Z Zy 


M 


Zt 


2,0,3 

Zq Z Zy 




^ p 

^(r(r^(rqf(r6-(r(r^(r(r^(D 

k,c 




qi,ai 


q2,a2 


^ P 

^ ,,, , ^ k,c 


qi,ai 



q2,a2 


Figure B.2.1: ^ 2 , 0,3 “direct” contributes to second order in opacity while ^20 3 

tributes to first order in opacity [10]. 


B.2 


,0,3 


The results that these calculations are based on appear in Appendix C in [10]. The 
relevant diagrams are shown in Figure B.2.1. The same procedure as described for the 
diagrams is used to determine that the matrix element for ^ 2 , 0,3 is given by 

=/‘ - 9 . 

X i l\M {p + k - qi - q2){-i) Amg {k - qi - q2){-i) Amg {k - 92 ) 

«J{p + [(c, a- 2 \,ai] (r„(2)r„(1)) x 

r d(^ r dq^ 

J 27r J 27r ((p + /c — gi — ^ 2 )^ — + ie) 

^ _ xg 2 1) 

{{k - qi- q2Y -mj + ie) {{k - ^ 2 )^ - + ie) ’ 

In order to evaluate the qf integral, it is convenient to rewrite the following 


g-*<?i(2:i-2^o)g-*<?|(2:2-2:o) _ g-i(gf+g|)(2:i-zo)g-*ij|(2:2-2:i) 


The first longitudinal integral (keeping with the numbering in [10]) 


h{p,k,qi,q2,zi - zq) = 


dqf v{Qi,qi)e ^ 0 ) 

27r {{p + k — qi — 52 )^ — + ie) 

1 


{{k -qi- 92 )^ - m 2 + ie) 


(B.2.2) 


In the same way that zi — zq > 0 determined that the contour be closed in the lower half 
of the qf plane, we again close the contour below the real axis here. The poles are very 
similar to the poles for and can be obtained quite simply by letting qi ^ qi + q 2 
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or making a substitution into the results of = qi + q 2 which gives the poles as follows 


and 


q^ = -Wo -ujm-ie 

^(1) . z ~ ■ 

(ll + 92 = -^0 -OOm-ie 

= -gf - Wo - Wm - ie, (B.2.3) 


= -Wo + W3 - ie 

2(2) I ^ I 

91 +92 = -Wo + W12 - ^e 

9 i^^^ = -92 - Wo + Wi2 - ie (B.2.4) 


The third pole, from the potential, does not need any such a treatment and therefore has 
no gl dependence. We therefore have that the poles are at: 


2(1) 

9i 

2 ( 2 ) 

9i 

2(3) 

9i 


= -q2 -uo-Com-ie 
= -ql - Wo + W12 - ie 
= -i/ii 


The residues are therefore also similar. From Section B.1.3 we have that 


Res [/2,9 i^^^] 


-v{-q^ - wo - Wm, q 

x{P+Y{Gjm + OJ12) 


Res 



x(P+) 2 (Wm + W12) 


(B.2.5) 


(B.2.6) 

(B.2.7) 


The residue due to the pole at = —i/ii is given by 


Res[gJ^*'^^] = lim 




( 3 ) 


u(9f, qi)((7f + i/ri)e ^ 0 ) 


{{p + k- qi - q2Y + M2 + ie) {{k - qi - 92)^ - + ie) 


(B.2.8) 


The individual parts of which are calculated in sections A.7.2, A.7.2 and A.7.3 so that 
the residue of this pole is given by 


Res[qf) = -m] = - 

—2ip,\ P+ q 2 — ip,i xP+ gf “ 


dvra, 1 


1 


1 


1 


dvra. 


(2z/q) x(P+)2 (g| - i^i)2 


„-i(g|-i/^l)(2i-2o) 


(B.2.9) 
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We can therefore calculate the integral I 2 (keeping in mind that the contour is once again 
negatively orientated) 


hip, k,qi,q 2 , “ ^ 0 ) = -^ 


1 


_x(P+)2(a;m + W 12 ) 

X (vi-ql -U0 + uu, q,)e*(-o-i2)(.i-.o) 
- vi-ql -uo- oj^, q^)e*(-o+-^)(.i-.o) 


dvra. 


1 


1 


(2i//i)x(P+)2(g|-i^l)2 
1 




x(P+)2 


U!m + W12 


{vi-q 2 - Wo -ojm,qi)e 


i{ujo+oj m ){zi-zq) 


- vi-ql - wo + wi2, qi)e*(‘^“-‘^i2)(^i-^o) 

dvra, 1 


+ 




2z^i (gf - i(Ui)2 

^(^2 d~ <^W, q;^) _ gi(a;o-i^l2)(zi-2o) 


x(P+)2 

dvra^ 


+ 


Wm + W 12 
1 


2z^i (g| - i(Ui)2 


-i((j|-iAti)(2i-2o) 


(B.2.10) 


Notice here that the potential is evaluated near which still needs to be integrated over. 
In light of the complexity of what follows, in order to be absolutely explicit, I will rewrite 
the Matrix Element here with the calculated h- So, the matrix element becomes: 


^ 2 , 0,3 =Jip + [[c, a 2 ], ai] (T,, (2)T,, (1)) x 


X (-*) 
X 4Euj 


d^qi 


H) 


rf^q2 02 


2igse • (k - qi - : 


+ 


x(P+)2 

dvra.. 


(27r)2 ' J (27r)2 

dQ2 u(g2, q2)e“*'^2(22-2o) 

(27r) (k - g' 2 )^ - m2 + ie ^ 

'^(^2 d~ <^W,q]^) _ ^i(uJo-iz!i2)(zi-zo) 


Cjm + W12 


+ 


2i/ii (g| - ipLiY 




(B.2.11) 
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A note about 5uj: We can ‘factorize’ the potentials out if we take the lo contributions 
to be the same. This has to do with the fact that, within the potential, w’s get neglected 
w.r.t. the /i (which is why the u(a;i,q]^) ~ u(0,q;^)). And so, for any pole that contains 
u, we have that 


v{-ql -oJo-Co, -uo-Co, 


u(0,qi)[e 




Similarly, if the pole contains fi, the /x dominates the ox’s, so we have 




^i{LLJo-coi 2 ){zi-zo) _ i(aJ 0 +w)( 21 - 20 ) 


The question might arise as to whether or not this is still valid for other poles (being the 
ones that will be kept at a later stage) and some reflection will reveal that it does indeed 
since those poles were suppressed exponentially due to {zi — zq) and not oj. Now we can 
more clearly see that the factorizing won’t happen as easily as before. So let us start by 
defining: 


W^e 


i(uJO+uj 

m. )(z\-zq) 


i(cjO-uJl2)(zi-zo) 
^ ? 


(B.2.12) 


which allows one to redefine I 3 . 


h{k,qi,ci 2 ,Z 2 - zi) 


f dq2 v{q2,ci2)e *'? 2(^2 2 : 1 ) 

J ( 27 r) ((/c - 92)^ -mj+ ie) 

'v{q2 +(iw,q;^){*} dTra^ 

d)m + wi2 2i/xi (gf-Vi)^ - 


(B.2.13) 


In the most general case (that includes the contact case), there are three contributing poles 
(two from the potentials and one from the propagator). The poles from the potential are 
calculated in section A. 6 . These poles contribute regardless of the system size, because 
the suppression in the large system approximation occurs due to the fact that {zi — zq) is 
very large. Here, this does not hold sway, since we will in fact take {z 2 — zi) to become very 
small. The terms in 5oj are neglected in the calculation. The pole from the propagator is 
calculated in section A.5.3. The only pole that is new in the small system limit is 
+i/xi which is not considered because it lies in the upper half of the complex 
plane. We therefore have the poles 


^ 2 ^^^ = —ifii + 6 u! (B.2.14) 

= -if^2 (B.2.15) 

^ 2 ^^^ = 0 J 2 — eoQ — ie (B.2.16) 

= +il^i, (ignored) (B.2.17) 
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Note however that, although there are no extra poles, there is an extra term, which will 
change the other poles. 

z(l) 

We calculate the residues individually: For ' = —ilJ-i there are three contributions to 

h-. 


Res = lim 


^ z(l) 


v{(lh q 2 )e 

((^ - + 'm?„ + ie) 


X 


a 


v{q2 + (^a;,qi){^} dvros e 

Wm + W12 2z^i (gf - Z/Xi)2 


(B.2.18) 


Since this pole comes from the potential inside the square brackets, the second term 
will evaluate to zero here because it does not contain anything to cancel. We therefore 
have almost exactly the same pole, we just include here {*}. Using sections A.7.3, A.7.2 
and A. 7.3 we have that The first residue is not affected by the small separation 
distance approximation: 

2 ( 1 ) _ 1 

L 3,92 ^p+ -2im {nl - fij) COm + ^^12 

(B.2.19) 


For ^2 


2 ( 2 ) 


—ifj .2 we have the following contributions 


Res [Is, ( 72 ^^^] = lim 


7 ■ 2 ( 2 ) 


^^(g2,q2)e ""^^92+^^ 2 ) 

({k — qY — rUg + ie) 


u(92 + ( 5 a;, qi){*} dvra., e */^i)(2i ^o) 


Urn + a;i 2 


+ 


2i/ii (g| - i7Xi)2 


(B.2.20) 


The results from sections A.7.3, A.7.2 and A.7.3 allow for the calculation of the second 
residue: 
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Res 


47ras 1 e ^^ 2(22 21 ) 

-——;-r— X 

-2i/U2 xP+ {-iM2) 

47ras r ^ 47rQ;^ 

^ I,,,, I, ,, I - - - __ T "A" f _ _ ^ „ 

.(^1 - + ‘^ 12 ) ^ 2 ifJ-i (^ 2 -/^ 1 )^ . 

(47ra^)^e“^2(z2-2i) 

(- 2 ^i)xP+ ^ 


{*} 

- /i 2 )(Wm + W 12 ) 


g-(M2+*Ml)(2l-2o) - 

( 2 ifJ,i)(fl 2 - /Ui )2 


(B.2.21) 


Here we see again the effect that cancels terms, causing the second term in square brackets 
to be suppressed as compared to the first (Section A. 10). Therefore, the second residue 
remains unchanged in the small system approximation. 


Res [h,qf ^] = 


(47ra^)^ {^} _ 

2^2 xp+ (^2 


g-/i 2 ( 22 - 2 l) 

- uDiUm + UJ 12 ) 


(B.2.22) 


2: ('3') 

The third residue is probably the simplest, so for qf' ’= 0 J 2 — ^0 — i^ 


D rr R3)i 1 . iv{q^, q2)e - a;2 + wq) 

Res 13,92 = \ - 771 -7^1-5-rx- 

1 [{k — 9)2 — m 2 + le) 


95-5-92 


u (92 + (fw, qi){*} ^ dvros e *•^'^2 zo) 


dim + <7712 


2 i^i ( 9 I - imf 


(B.2.23) 


where we have, using the same sections as for the previous poles that the residue due to 
the third pole (from the propagator) is then 


dvra^e ^^o)(22 21 ) 

Res [73,92' ' = W2 - WoJ = - X 

dvra, 


xP+ 

dvra^l*} 


.(/ii)2(a;m + W 12 ) 2iiJLi 
{4TTasf e-*(‘^2-^o)(22-2i) 


{lA) 


xP+ 




{*} 


Ulr] 


WI 2 


g —ZO) 

2 i//i 


(B.2.24) 


We see again the suppression effect from Section A. 10, suppressing the second term relative 
to the first. Therefore, the third residue remains unchanged under the small separation 
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distance limit. 


Res = 


(dvrag)^ e ‘^o )(^2 zi) 


xP+ 




Wm + Wi 2 


(B. 2 . 25 ) 


We now need to sum these terms in order to calculate the integral I3. 

1 47 ^ 0 ^ Attus r , 

=Spr—i^r—MW?) {*} 


(dTra^)^ {:*:} 


0-/22(22—21) 


+ 


2/r| xP+(/if - ^|)(a;m + W 12 ) 

(dvras)^ g-*(‘*^2-t2>o)(z2-2i) 1^1 


xP+ 

(47ra2)^(-^){*} 
xP~^ (^LOrn “ 1 “ ^ 12 ) 

1 

+ 


lJxl4 


+ ^^12 
g-i(aJ2-i2;o)(22-2l) 

/^ 1^2 


2(/^2 - /^l) 


r g-Mi ( 22 -^ 1 ) 

g-/12(22-2l) 'I ■ 

1 Pi 

^2 i. 


(B. 2 . 26 ) 


This seems difficult, but we need only consider two extreme limits (ignoring for the mo¬ 
ment a common factor of - —): 

1 . The limit of well-separated scattering centres Z 2 — zi 1 /p.] 

2 . The special “contact” limit 2:2 = zi limit to compute unitary contributions 

For the first case in which the scattering centres are separated on the scale of the mean 
free path (z2 — 21 ~ A S> 1 //^), equation (B. 2 . 26 ) reduces to 

h{k, qi, q 2 , 22 - > l/p) qi)^^(0, q 2 )e“*(“ 2 -^o)( 22 - 2 i)^ (B.2.27) 

where the other terms have fallen away because the exponential terms have massively 
negative arguments. For the second case in which Z 2 — zi = 0, consider the following fact 
that arises from pf = pi2_ = qf -|- /i^ 


{4TTasf / 1_1 \ _ _ 


(dvras)" 


2(i^2 -/^i) V/^2 lAJ 2(/r2 - qf -/i2 q2) p"^ - - qi 


1 


1 


(dvra* 


2(qi - q?) 


(qi - q?) 


(/r 2 -q 2 )(/r 2 _qi) 


= 2^(o>qi)^(o,q2)- 


(B.2.28) 
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We therefore have that equation (B.2.26) reduces to 


i 1 


h{k,qi,q2,Z2 = zi) ^ qi)u(0, qa) 


(B.2.29) 

Which is exactly half (in strength) of the first case. 

We can now return to the amplitude for ^ 2 , 0,3 in the contact limit. Note first that 


AEoj 


1 


4Eoj 


1 


!:(P+)2(a;(i2) + uj'm} 2k+ P+ [(k - qi - q2)^ + + m^] 2/c+ 

^_1_ 

[(k - qi - q2)^ + + m^] 2 ^;kE^ 


1 


1 


2 [(k — qi — q2)^ + + m^] 


We therefore have the amplitude: 


.^ 2 . 0,3 s3J(p)e*'''+''>«(-i) j j ^„{ 0 . x 


d?‘q2 


(27r)2 

X ^ 2^2 + ^2] «2], ai] {Ta,T^,)x 

X |gd“0+<^m)W-2o) _ gi(tJo-‘^(12))(^l-2o)| 


(B.2.30) 
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M 


2,0,0 

Zo 

^ • i 




k,c 


Z-) 


% 




I 

® — 

q 2.^2 


2,0,0 

Zl 


k,c 



Figure B. 3 . 1 ; ^200 graphs for the well separated case as well as the contact limit 

[ 10 ]. 


B .3 ^2,0,0 and ^2,2,0 

B.3.1 .^ 2 , 0,0 


The Matrix element for .^ 2 , 0,0 is calculated as before from the diagram in Figure B.3.1 

^ 2 , 0 , 0 =/ (2"J)4 + fc - gi - X 

V{q2)e^^^^H-i‘^Efi9s{2p + k)^e>^x 

X i Am {p + k - qi - q2)i Am {p - Qi - 92)1 Am {p - q2)aic{Ta^Ta^) 


j(p)edp+fe)^o 

H)^ 


d\2 „-iq2-b2^ 
(271)2 X 

X 

2igs{e 

X 

•k) 

a2aic{Ta2Tai){2E)‘^ X 



f dql 

dq2 





J 27r 

27r 

((p4 

-k - qi - 52 )^ - + 






t)(q^)e~*'^2(^2—^o) 



{{P- 

91 - 

■ 92 )^ 

- M2 + ie)((p - 52 )^ - 

- M2 + ie) 


Here we define (again keeping with the numbering in [10]) 


l 2 (p,fc,qi,q 2 ,zi - zo) 

^ fdqf _ u(gf,qi)e^(^i-? 2 )(^i-^o) _ 

J 2tt (^{p + k - qi - q 2 )‘^ - + ie) {fp - qi - ^ 2 )^ - + ie) 


There are three poles, one from each propagator and one from the potential. These are 
found by (as for .^# 2 , 0 , 3 ) making the substitution gf = gf + i?!) using the results from 
sections A.5.2 and A.5.1 and then ‘solving’ for qf. This results in the following three 
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poles: 


= -g| -uo-Com-ie (B.3.3) 

= -Q2 + x{uj(^i2) - Wo) - ie ~ -gf - ie (B.3.4) 

= -ipii (B.3.5) 

zil) 

We calculate the residues individually. The residue due to the pole iAqi=qf' is given 
by 


Res [12,91 
= lim 


ti) 


■1^(91, ^o)(g2 _g _g _g (2;^) 

((p + fc - 9i - 92)2 - M2 + ie) {{p - 9i - 92)^ - M2 + ie) 


(B.3.6) 


Sections A.6, A.7.2 and A.7.1 give that the first residue is 


Res [/2,9 i 


'^Til ~ '^( 921^1) i((^n+d!rr,.)(zT-Z()) 


(p+)2(_^0 - Wm) 


(B.3.7) 


The second residue due, to the pole at 91 = —9! + x(a;(i2) — wq), is calculated as follows: 
Res [/2,9 i^^^] 


= lim 

7 2(1) 


u(9f,qi)e d^f+gfKzi 2 o)(-gj_^g| 

((p + fc - 9i - 92)2 - M2 + ie) ((p - 9i - 92)^ - M2 + ie) 


(B.3.8) 


The individual contributions are calculated in sections A.6, A.7.2 and A.7.1 to give 


Res \h ~ 1 Q 2 ^^i) i(o)(zi-z2) 

Res[/2,9i + 

1 u(-9f,qi) 


(p +)2 _|_ 2 )^^ 


(B.3.9) 


The residue due to the pole at q^^^^ = —ipi is given by 

Res [/2,9 i^^^] 

= lim 


91^91 


( 3 ) 


^( 9 i,qi)e 2 :o)('gZ | 


((p + /c - 9i - 92)2 - M2 + ie) ((p - 9i - 92)^ - M2 + ie) 
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with simplifications from the same sections giving 


Res 


-¥i] 


- 2^1 (P +)2 - i^i )2 


(B.3.10) 


We can now calculate the integral l 2 - 


h{p,k, “ ^o) 


—2'Ki 

27r 


At: a. 


ai{u)o+u) 

m ){zi-zo) 


+ 


+ ^l) (^+)^(-W0 - OJm) 
Alias 1 1 

+ ^l) (‘^O + djm) 


Anas 1 

- 2 i/ii (P +)2 ((^1 - 


(B.3.11) 


We also see here the suppression effect, so the second term is again suppressed relative to 
the first, and the integral stays the same: 


/ 2 (p, A:,qi,q 2 , 2 ;i - ^o) 


(47ra^)(i) _ i 

(P +)2 (wo + Wm)(g '2 - 


(B.3.12) 


The matrix element with the calculated l 2 - 

^ 2 , 0,0 J ^e-^^^-^^{-i) J 

^'^^^^'^K 2aic{Ta,Ta,){2Efx 

dq2 u(q'"2)e“*''2(^2-^i) (47rQ;s)(i) - 1 

^ 27 : {{p - q2)‘^ - M'^ + ie) (-P+)^ {utq + dJm){q2 - 


The remaining integral, I 3 , is then defined as 


f dq^v{-q^,q^)v{-q^,q2)e 

*(P. q., q,. - . 0 ) = / ^- (fa M2 - 


{{p - 92)^ - + ie) 

^i{aJO+u> m ){zi-zo) _ I 


(wo + dJm){q2 -iPlY 


(B.3.14) 


Since we will specifically be looking at the contact limit, the poles from the potentials can 
not be neglected. We therefore have three poles: two from the two potentials and one 
from the propagator. Also, there are no new poles from the short path length limit, and 
therefore the rest of the calculation of this diagram is identical to the large path length 
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calculation. By sections A.6 and A.5.1, these poles are 


^( 1 ) 

Q2 = -*/^l 

z(2) 

92 = -*/^2 

z(3) (k + q)^ - 
92^ = -- 


(B.3.15) 

(B.3.16) 

(B.3.17) 


z(S) 

For the last pole, it will be sufficient to use q 2 ^ ' ~ —ie, but the full results here for 
completeness until the residue for that pole is calculated. Again, we calculate the residues 

z(l) 

individually, starting with the residue due to the pole at ^ 2 ^ ^ given by 


Res [13,92^^^] 


lim 


^^(-92^qi)^(-92>q2)(92 +^m) 

{{p - 92 )^ - + ie) 


and computed using sections A.7.2, A.7.3 and A.7.1: 


(B.3.18) 


Res [13,92^^^ 


-ipi] 


{A-Kas) {AttUs) i 
-2i^i - p,l P+m 

-(47ra,)2 

2 /r 2 (;. 2 _^ 2 )p+ 


(B.3.19) 


The residue for the pole at = —ip 2 is calculated in much the same way with minor 
changes. It is given by 


Res [13,92^^^] 


lim 

l2^l2 


( 2 ) 


^^(-92^qi)^^(-92^q2)(92 +^M2) 

{{p - 92 )^ - + ie) 


and its parts are simplified with results from the same sections to give 


(B.3.20) 


Res [13,92^^^ = 


-ip2\ 


-(47ra,)2 

2 ^iP+(^ 2 _^ 2 ) 


(B.3.21) 


The last pole is at q^^^^ = —- — ie and gives the residue 

r^( -92,qi)^(-92,q2)(92 - 

{{p - 92 )^ - + ie) 


Res 

T -^(3) 

h,q 2 

= hm 



z 2(3) 

q^^qf L 


(B.3.22) 


the individual contributions to which are also from the above section and give 
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Res 



(k + q)2-kl 
P+ 


{Airas) ix{ui2-uj())(z2-z{) 


(B.3.23) 


Which means that the integral I 3 becomes 


h{p,k,q;^,q 2 ,zi - zq) = 


-(47ra,f 

^z;(0,qi)^;(0,q2)e--(-^-°)(^^-^i) 


—i 

P+ 

+ 


r;( 0 , qi)z;( 0 , 


(47ra^)2 / 

/g-Atl(z 2 -Zl) 

g-lt2(z2 —21) \ - 


V Ml 

M2 /- 


(B.3.24) 


Again we consider two cases, the well separated and the contact limit. 
1. Z 2 — zi ^ l/fi: The two terms in the round brackets go to zero: 


—i 




(B.3.25) 


2. Z 2 = zi: This case is slightly more complicated 


/3 = 


—I 

P+ 


u( 0 ,qi)u( 0 ,q 2 ) + 


(dvra^)^ f 'i- _ 


2(1^1-/^i) VMi PlJi 


(B.3.26) 


We can see that we get the factor of 1 /2 again in the following way. Consider 


(dvrcTs)^ / I 1 \ _ (4Tras)‘^ / I 1 \ 

2(^2 _ ^2) (^^2 ^ 2 ; - 2(q2 _ q2) [^2 _ q2 ^,2 _ q^j 

^ (47ra,)2 / \ 

2ie^^^{{n^-ql){p^-ql)) 

= - 2 ^(o>qi)^'(o,q 2 ) 


(B.3.27) 
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We therefore have that 

h ~ (B.3.28) 

Finally, the amplitude can now be written by combining the above results (remembering 
that part of the result for I 2 is integrated over in I 3 and so should not be included twice): 


^ 2 , 0,0 f 


X ?^^^^a2aic(r,,r„J(2E)2x 


dqf dq. 


(27r)2^ ^ 'J (27r)2 

,^Ta,){2Efx 


J 27r 27r {{p + k — qi — 92 )^ — Af 2 + ie) 

u(q-^2)e“*''2(^2-^o) 

{{p - qi- q 2 Y - Af2 + ie) ((p - q 2 Y - M2 + ie) 

d^qi r d^q2 _ 


i{p+k)xo 


(27r)2 J (27r)2 


ij{p)e 

X • k)a 2 aic(raira 2 )( 2 £^)^x 

X 

X ^ _ ^ f i(uJo+i:jm)(zi-zo) _ x'i X 

(P+)M^o + UJm) ^ 

(—i) 1 , separated 

X ^^t’(0,qi)?^(0,q2) X <j ^ 


*(qi+qi)'bix, (bi Ribs) 


2 ’ 


contact 




X 2iff^(e-k) ^^i(a;o+i:,^)(2:i-zo) _ ^ 


k^ + m2 + M2 x2 


X a^axcTd^Tai 


1 , separated 
^, contact 


(B.3.29) 


(B.3.30) 

(B.3.31) 


(B.3.32) 

(B.3.33) 


Where I have used that 


1 


x( 2^)2 

x{p+Y 

1 


= 1 


and 

1 


xP+ CUo + dJn 


k^ + m2 + M^x^ 


(B.3.34) 


Therefore, this diagram, ^ 2 , 0,0 remains unchanged under the small separation distance 
limit. 
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2,2,0 

Zo 


Zl 
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qbai 


k,c 








I ' 

q 2, a 2 



Figure B.3.2; ^ 2 , 2.0 diagrams in both the well separated and contact limit case [10]. 


B.3.2 .^2,2,0 

Figure B.3.2 shows the diagram for ^ 2 , 2 , 0 , the amplitude for which is very similar to 
that for ./#2,o,o and differs only in the propagators (and the order of the color generator 
matrices in the result): 




(27r)4 (27r)^ 

X x {-i2Efigsi2p + k)f,e^^x 

X i Am (p + k - qi - q 2 )i Am (p + k- q2)i Am {p + k)aic{Ta2Tai) 


p,), 

X ca2ai{Ta^Ta^){2E)‘^x 


(271) 


dqf dq. 


u(q^)e 


J 27t 27r {{p + k — qi — ( 72 )^ — + *e) 

u(q^)e~*'^2(^2-^o) 

{{p + k — q 2 Y — + ie) {{p + k^ — + ie) 


(B.3.35) 


Notice that the last propagator does not depend on qi or 52 which greatly simplifies the 
calculations below but must be kept in mind in the final calculation of the amplitude. We 
define the first integral 


22(p, A:,qi,q2,^i - 2 : 0 ) 


dqf v{qi,e{i)e d<?l+9|)(2:i 2 ^ 0 ) 

277 ((p + A: — gi — ( 72 )^ — + *e) 


(B.3.36) 


There are two poles, one from the propagator and one from the potential (that had 
previously been neglected). Using the results in section A.5.2 along with an appropriate 
substitution in much the same way as was done to obtain equation (B.2.4), one can 
calculate the pole in the propagator, while section A.6 describes the calculation of the 
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pole at the potential. The poles in question are then at 


z(l) z 

Ql = -<l2 - Wo - Wm - ze 
z( 2 ) 

Qi = -*/^i 


(B.3.37) 

(B.3.38) 


The residue here is very similar to the corresponding residue in the amplitude of diagram 
but does differ slightly: 


Res = -gf - Wo - Wm] 

v{qf, (2)^^ 


= lim 


( 1 ) 


{{p + k- qi- q2Y - M 2 + ie) 


(B.3.39) 


The individual contributions are calculated in sections A.6 and A.7.2, giving 


Res [h, qf'^ = -ql - wq - w^] = 


(B.3.40) 


The residue due to the second pole at = —ifii is given by 


Res [l 2 ,ql^‘^^ = —ifJ-i] = lim 


i z(2) 


{{p + k -qi- q2Y - m 2 + ie) 


(B.3.41) 


We can again calculate each part individually using results from the same sections, re¬ 
sulting in 


^ rr R 2 )n dvra, 

Res / 2 ,gi = - 5 x 7 ^-^^ 

-2z7i P+(g|-z/ii) 


(B.3.42) 


Adding these two poles, the integral I 2 is given by (with the negative sign again because 
the contour is negatively orientated) 


d2(p,/c,qi,q2,2i - zq) 


2'Ki 


dvra gi((Jo+<^m)W-zo) 4 ^q, ^-i{q^-i^l+{z^-zo) 


.(g|)2 + 7f P+ 2i7i P+(gf-i7i) 

—i dvra rgi(ajo+Wm)( 2 i-^o) g-i((j|-i/ii)(2i-2o)- 


P+ (g| - iixi) 


gf -h iiii 


2ifii 


(B.3.43) 
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The matrix element with the calculated I 2 is given by 

f „-^qo.b,2ig.(e-k) 


, 2,0 ^Jip)e^^^%-iy 


(2vr) 


^ (2vr)2 
X ca2ai{Ta2Tay{2Ef X 

f dql 


J 27r ((p + /c — g' 2 )^ — + ie) {{p + + ze) 

—i dvra rg*(ajo+wm)( 2 i- 2 o) g-dQl-^z^iK^i-zo) 


P+ (g| - i/zi) 

So we define the following I 3 : 

/3(p,A:,qi,q2,zi - zq) 


gf + z/ii 


2ipLi 


f dql 1 

■gi(tJO+‘^m)(2l-2:o) 

^-i{q^-ifMl){zi-zo)- 

J 2-71 (gf - Z^l) 

92 + 'i-Pl 

2z/zi 

A-iras 

g-ig|(22-2i) 



(g |)2 + ^2 ^ _ ^^^2 _ m 2 + ze) 


(B.3.44) 


(B.3.45) 


Here the poles from the potentials can never be neglected because their contributions are 
not even exponentially suppressed in the large path length limit. We therefore have three 
poles, two from the potentials (see section A. 6 ) and one from the propagator (section 
A.5. 2). (The new pole has a positive sign and is therefore ignored because it is in the 
upper half of the complex plane. However, there is a new term, and so all the poles have 
to be recalculated. The poles are therefore 


z(i) 

92 = 

z{2) 

92 = 

■ 2 ( 3 ) 

92 = 


-IpLl 

-ip2 


(B.3.46) 

(B.3.47) 

(B.3.48) 


We calculate the residues in the same way as before - individually and by considering each 

z{l) 

contribution. The residue due to the pole ai qf' — ipi is calculated as 


Res [ 13 , 92 ^^^ = 


= lim 

2 2 ( 1 ) 

47ra, 


f ( 9 I +W) 

■gi(tJO+<^m)(2l-2o) 

g-dg|-*Mi)(2:i-2o)- 

1(92 - W) 

92 + *7^1 

2z/ii 


3 -i(}|( 22 - 2 l) 


( 92 )^ + pliijp + k- 92 )^ - M2 + ze) 


X 


(B.3.49) 
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Each part can be evaluated individually using equation (A.7.3) and section A.7.2, resulting 
in the residue 

2(1) _ . ^ _ (dvra^) ^i{a:o+u^)izi-zo) 


Res [73,^2^ - -i/Ui] 


^2 - lA -2ifii 

m ){zi-zo) 

2 P+(/ri-/r2)^2 


(B.3.50) 


This is exactly the result for the large path length calculation, the residue due to the pole 

, 2 ( 2 ) ... , 
at 52 = IS given by 


Res [73,92^^^ = -*/^ 2 ] 
= lim 
dvra^ 


f (^2 + *(^ 2 ) 

■ ^i(u)o+^m)izi-Zo) 

^-i(q2-iU\)Dl-zo)- 

l(5f - im) 

gf + ipi 

2ifii 


,-iq^{z2-zi) 


X 


{Q 2 A + IJ'liiP + k- 52 )^ - M 2 + ie) 
Which can be calculated using the same results as above: 


(B.3.51) 


Res [ 73 , 52 ^^^ = -iP 2 \ = 


Ana,, e zi) 


(-2i^2) -i{pi + P2) R+(-i/i2) 

gi(a;o+iI>m)(2i-2o) g-(Mi+M2)(2i-^o)' 


i{pi - P 2 ) 


—2ifj,i 

,-{lll+fl2){zi-Zo) i{uJo+i:>m)(zi-Zo} 


2P+nl{ni + ( 12 ) [ 2^1 (/ii - ij,2) 

z( 3 ) 

The residue due to the third pole at 52 ^ ^ = “^0 ~ — i^ is given by 


(B.3.52) 


Res [ 73 , 52 ^^^ = “^0 - Wm] 


= lim 


f (^2 + ^0 + ^m) 

■gi(a;o+iI>m)(zi-2o) 

g-i(i}|-i/ 2 l)( 2 :i-zo) ■ 

1 (^2 - */^l) 

g| + i/xi 

2i/Ui 


dvra, e-*'' 2 (^ 2 -^i) 


( 5 f )2 + + 52 )^ - M 2 + ie) 


X 


(B.3.53) 
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with the following contributions from sections A.7.2 and A.7.2 to give 


z(3) A-KOis 1 ^ ^i{ujQ+Gjrn)(z-L-Zo) Zq)- 

fil P+ -ip,i 


Res [13,92 ] = 2 

h 


ip .1 


2 ipLi 




^i{bJO+Cjm){zi-Zo) _ bg-Ml(21-20) 


(B.3.54) 


We combine the residues (keeping in mind that the contour is negatively orientated) to 
obtain the result for the second integral I 3 . 


2'Ki f (47ras)e zi)gi(ujo+u!m)(zi zq) 

/3(p.t.q,.q..*, - 3o) = -^1 -- 


+ 


47ra e“^2(2;2-^i) rg-(ui+U2)(zi-zo) gi(ajo+i^m)(2i-2o) 


+ 


2 P+^ 2 ( 1^1 + 1 ^ 2 ) 

47rQ;oe*(^°'’'‘^ m )(22-Zl) 


2 /ri 


(^1 - M 2 ) 


i(uJo+cjm)(zi-Zo) _ ^ -^ 11 ( 21 - 20 ) 

2 


-P+M 1 M 2 

(-i)47ra, f e-m(22-2i)gi(a;o+i; m )(2l-20) 

i 2 (m 1 -m?)m? 

q-U2{z2-Zi) / ^-{ni+fJ,2){zi-Zo) gi(^0+‘^m)(2l-2o) 


+ 


2 M 2 (Mi +M 2 ) 

gi([2o+iI> m )(22-2l) 

M 1 M 2 


2 mi 


(M2 - Ml) 

1 


gi((3JO+<l’m)(2l-2o) _ _g-Atl (21-20) 


(B.3.55) 


Before we can really look at the contact and large separation cases, it might be pertinent 
to simplify the curly brackets. So, let us define: 


g-Ml (Z2 -2i) gi(t,JO+“m) (21 -20) 

■ 2(^i-Mf)Mf ■ 

g-/12(22-2l) / g-{fil+U 2 ){zi-Zo) gi(t^O+‘l’m)(21-2o) 

2/r|(;Ui +/X 2 ) V 2/ri (mi - M 2 ) 

i(‘^0+Wm)(22-2l) ,1 , , 

_|_ I_ I ^IpO+l^mjlzi-Zo) _ ^^-/ll{zi-Zo) 

M 1 M 2 V 2 


(B.3.56) 
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Which would give the matrix element 




^.-*qrbi 

( 2 vr)^ 

X ca2ai(ra2raj(2^) 

=J{p)e^P^^{-i) 


(27 
dvra^ 


(-7)47ras 


(27r 


X ca2ai(ra2rai)(2^)^^ 


-P"*" ((p + /c — ( 72 )^ — + ie) P~^ 

^.-iqrbl/_7^ f -iq2.b, 2ig.(6-k) 

^ V (27r)2® x 


(P+)2 777,2 _|_ 15-2 _|_ ^2]yp2 


(-!)■ 


(B.3.57) 


To simplify this, let’s look at {<8)}: 


(-i)M=< 


j-/il( 22 -^l)gi(‘^ 0 +‘^m)( 2 l- 20 ) 


_p-M2(22-2l) 


2 M 2 (P 1 +P 2 )L (7ti-li2) 

gi(a;o+(2' m )(22-2l) 


^i{uiO+i^m)izi-Zo) g-(/il+/i2)(2l-2o) 
^ 2//1 




1 


^i(oJO+Ulm)(zi-Zo) _ _g-7tl(zi-2o) 


3-/J.I (22-21 )gi(LJ0+<2)m)(21-20) 


+ 


q-U2{z2-Zi) ^i{u}Q+Cjm){zi-ZQ) g-/12(22-2l)g-(/ll+/12)(2l-2o) 


271 ^( 7 x 1 + 7 ^ 2 ) (M 1 - 1 ^ 2 ) 2pLl{pLl + pL2)2pii 

gi(a;o+(im)(22-2l)gi(lJo+<2'm)(2l-2o) 2 g4‘^0+‘2’m)(22-2l)g-/ll(2l-2o) 




+ 


i4i4 


(B.3.58) 


Here there are two terms that appear in the small separation limit that did not contribute 
to the large separation distance calculation. The last term on the second line and the last 
term on the third line are both additional contributions that are exponentially suppressed 
in the large path length limit. From this we also see that this expression still reduces to 
the large path length result. One interesting thing that arose in the large path length limit 
was a factor of 1/2 that appears in the contact limit. We can check that this factor still 
appears if we take this expression for the short path length and look at the contact limit 
as well as the large path length limit. The two ‘new’ terms will be zero (large path length) 
while the exponents containing {z 2 — zi) —)■ (zi — zi) = 0 become 1 and the exponents 
containing (cuq + — zq) are factorized out, so that we have that 


(^_X)|(g) I contact ^gi(a;o+‘^m)(2i-2o) 


1 


+ 


1 


+ 


g-(All+A12)(21-2o) g-/ll (21-20) 

+ 


+ 7 x 2 ) 




1 ■ 


(B.3.59) 
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Now, the square brackets in the first line are exactly the terms that appear in the large 


path length limit and they still produce the factor of 1/2: 


1 1 


1 

l4l4 


= _ ^ - 

1 1 


_3_ 

1 ^ 2 / ^ 1^2 


1 


21 ^ 1^2 


(B.3.60) 


And so, although we are not surprised, we are relieved. Suppose we now only consider 
the expression for the short path length in the contact limit, then we have that 


(_l){(g,}contact ^ _ 1 _Lgi(a;o+<:- m ){zi-zo) 
2 Ml M 2 


Ml M 2 

g-Atl(zi-zo) 
-I-- 


g-(m+/i 2 )( 2 i-^o) 


(B.3.61) 


Now, we would like for these new terms to mean something, so notice that we have 

1 ^ 1 Ml 

4MiMi(Mi + M2) MiM2 4(mi + M2) 


We also have that 


g-(Ml+M2)(2l-2o) _ g-Ml(2;i-^o)g-M2(2l-2o) 


And so, the two small separation distance terms become: 


1 g-/il(2l-20) g-(Ml+M2)(2l-20) 

2 2|u|(/ii +/12)2 mi 


1 

2 Ml M 2 


^Ig-M2(2l-20)X 

2(mi + M2) / 
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M 
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Figure B.4.1; Diagrams ^ 2 , 0,1 and ,^# 2 , 0,2 showing only the special contact case [10]. 


Which finally gives the matrix amplitude for ^ 2 , 2 , 0 - 


J J X 


X ca2ai{Ta2Ta^){2Efx 


(-l)(47rag )^X 

- gi{uJQ+i:jm)izi-Zo) g-Atl(2l-2o) / 

filfll Hjfll \ 


^^g-M2(21-^0) Y 

2(^1 + /i 2 ) / . 


X 


{ 1 , separated 
^, contact 


(B.3.62) 


It is worth the effort to rewrite this expression to more closely match the form of the 
other two scattering centre diagrams since it will greatly minimize the work necessary 
when squaring and summing in the final energy loss calculation: 


^ 2 , 2,0 j 

X C(l2Ch\{Ta\ 


qi)ao, qi): 

, (- 1 ) 0 -k) „ 




J{uJo+i:)m)izi-Zo) _ ^-fll(zi-Zo) j I _ 


file 


-U2(zi-zo) 


2(/ii -I- p,2) 


{ 1 , separated 
I, contact 


(B.3.63) 
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B .4 ^2,0,1 

and ^2,0,1 

B. 4.1 ^2,Q,l 



Figure B.4.1 shows the diagram for the contact case for ^ 2 , 0 , 1 - The matrix element is 


r d^ d‘^ 

J (2^)4 + k-qi- x 


(27r)4 (27 

X (-zF;+)Aii Am {{p + k-qi- q 2 )i-i) Amg {k - qi)i Am {p - ^2) 

„-*qvbi f Mag- 


(2.)- y (2.)^ 

X 2i5,(e- (k-qi)e“a2[c,ai](r,2r„J(2S)2x 

f dqf dq 2 

J 277 27 r {{p + k — qi — ^2)^ — + ie) 

u(q' 2 )e“*'^ 2 (^ 2 -^o) 

{{k - qi)^ - m2 + ie) {{p - 52)^ - + ie) 




(B.4.1) 


We define the qf integral (keeping with the numbering in [10] as 


/ 2 (p, A;,qi,q 2 , 2 :i - 2 : 0 ) 

_ fdqf 


v{qi,(li)e *( 9 i +'? 2 )(^i ^ 0 ) 


27r {{p + k — qi — ( 72 )^ — + ie) {{k — qi)'^ — m^ + ie) 


(B.4.2) 


The three poles are calculated in sections A. 6, A.5.2 and A.5.3 (along with an appropriate 
substitution as was done to obtain equation (B.2.4)). These three poles are then at 


2 :( 1 ) 2 

qi = -q2 - Wo - Wm - 76 

2(2) , 

qi = -Wo + wi - ze 

2(3) 

Qi = 


(B.4.3) 

(B.4.4) 

(B.4.5) 


2 ( 1 ) 

The residue due to the first pole at ^ = —gf — wo — Wm — ie is calculated as follows 
Res [/ 2 ,gi^^^ = - Wo - Wm - ze] 

(B.4.6) 


= lim 


( 1 ) 


{ql + q 2 + OJQ + ujm)v{ql,eii)e *(9f+9|)(2i 20 ) 


{{p + k- qi - ^ 2 )^ - Af2 + ie) {{k - gi)2 - m^ + ie) 


the contributions of which are calculated in sections A.7.3, A.7.2 and A.7.3 to give 
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Res [l 2 ,ql^^^] 


v( q2 cup -i(ujQ+i:j^)(zi-zo) 

x{P+)‘^{q^ + uJi - Wo) 


(B.4.7) 


z( 2 ) 

The residue due to the pole at = — wq + wi — ie is calculated with 


Res [l 2 ,ql^‘^^ = -Wo Twi] 

{qf + Wo — wi)u(—Wo + wi, 


= lim 


{{p + k - qi - q 2 Y - M2 + ie) {{k - qiY - + ie) 


(B.4.8) 


Again we can most easily calculate this by looking at the individual contributions calcu¬ 
lated in the above mentioned sections to give 


Res [l 2 ,ql^‘^^ 


-Wo -h Wi] 


1 v{ Wo-|-Wl,q;^) 

xP+ p+{ijji 


(B.4.9) 


The residue due to the third pole 


Res = -ipii] 


= hm 

z , 2(3) 


{qf + iHi)v{ql,e{i)e 20 ) 

{{p + k- qi- q 2 Y - M2 -|- ie) {{k - gi)2 - -|- ie) 


dvra^ e“*('^M*A‘i)(2i-2o) ^-jrag e“*(^M*Mi)(2i-2o) 
-2ifii x(P+)2(-i/ii)2 “ 2 inl x(P+)2 


(B.4.10) 


We therefore have that the integral I 2 become 


^ 2 ( 4 ', A:,qi,q2,2:i - ^o) 


r u(-gf - cup - Om, qi) 

2tt x(P+)2(g|-I-wi - Wo) 
_1 u(-wo-4 wi,qi) 

xP+ p+(u;i -4 Wm + gf) 
dvra, e-*(9f-*Mi)(2i-2o)' 


-h 


2 ip,l x(P+)" 


{-^) 


v{ q2 Wq U)rl^,^.l) -i(uj(,+Lj^)(zi-Z(t) 


x(P +)2 [ (gf-4Wi-Wo) 

_ + ^ 1 ; Tl) ^i(q^+uJo+LOi)(zi-zo) 

(wi -|- Wm + <?!) 


-h 


dvra 


2 ifi'l 


£p-d9|-*Mi)(2i-2o) 


(B.4.11) 
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We also see here the major cancellation effect here, suppressing the small separation 
distance correction term in relation to the large path length terms. 


h{p,k,cii,q2,zi - zo) 


(-^) - ^0 - tT'm, qi) 

x(P +)2 (gf+cji-wo) 


v{ <-^0 + ^i(q^+CLJo+LLJi)(zi-zo) 

(wi + Urn + gf) 


(B.4.12) 


We now turn our attention to the integral which we define as 


h{p,k ,qi,q2,2:i - Zq) = J 


dq^ 


v{q2,^2) 


27r wi + Wm + gf - {{p - g 2 )^ - + ie) 

_ e-ff9f(^2-^o)-(^i+^o)(^i-2o)^(_^p + wi, qi) 

f dq^2 v{qh^2) 


27 r Wi + Wm + gf - ((p - g'2)^ - + ie) 

- cuo - cD™, qi) 




(B.4.13) 


Notice here the difference between the first term in the exponential in the second line and 
that of the first line. It is important to note that one reads {z 2 — zi) while the other is 
(zi — Zo)- Take care also to note that in the second exponential, the term {z 2 — zq) appears 
which causes an exponential suppression in the large separation distance limit. There are 
still only three poles here; two from potentials (note that the very last potential does not 
depend on gf and therefore does not contribute a pole), and one from the propagator. 
The singularity that arises in the first fraction is in the upper half of the complex plane 
and is therefore not considered. These poles are calculated in sections A.5.1 and A.6 to 
be 


z(l) 

Q2 = 

z{2) 

g 2 = -'>'P2 

z{3) 

^2 = -Wo - W„ 




-ipLi 


(B.4.14) 

(B.4.15) 

(B.4.16) 
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The residue due to the first pole at 


Res [Is, 92 ^^^ = -ie] = lim 


2 ( 1 ) 

Q 2 = 


v{q2,fl2)iQ2 + *e) 


p-i(j|( 22 - 2 l) 


uJi + Cjm + q2- {{p - q2Y - M 2 + ie) 

_ g-i(5|-^i+a;o)(2i-2o)^(_^^ ^ 


(B.4.17) 


The individual contributions of which are detailed in sections A.7.3 and A.7.1, giving that 
the first residue for this integral evaluates to 


Res [I 3 , 


2 ( 1 ) 

.92 


-ie 


(B.4.18) 


The residue due to the second pole (at 



—ip 2 ) is given by 


Res [ 13 , 92 ^^^ 


-ip2] 


lim 

2 2(1) 


^(92.q2)(g|+ ^A^2) e 

U}i+u}m + q2- {{p - 92)^ - M2 -h ie) 


X 


e*(‘"o+“-)(^i-^o)^(_g| - cuo - cD™, qi) 


(Anr/y p ^ ^ . 

_ V ^_ (^^{(DO-\-UJm)iZl-Zo) _ -fl2{zi-ZQ)\ 


(B.4.19) 


Where we have used the fact that too <€. P 2 - 
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z(3) 

The last residue is due to a pole at ^ 2 ^ ' = —wo ~ — ifJ'i- 


Res [I 3 ,q 2 ^^^ = -Wo - Cbm- «/Ul] 

^^(92>q2)(^2 + w) 


= lim 






Wi + Wm + gf “ ((p - ^ 2 )^ - M 2 + ie) 

i(a;o+< 2 >™)(^l-^o)^(_^| _ q^) 


dvra, g-uiiz 2 -zi) r 


mI -^1 


dvra.. 


— 2 ip,i 


^^i(uJO+d)m)(zi-Zo) _ _ffffs_^-tll{zi-Zo) 


— 2 z/Ui 


(dvra^)^ e ^i) 

Ml - M 2 ^+ 2 ^f 


gi(aJo+(I'm)(2l-^o) _ g-/il(2l-2o) 


(B.d.20) 


We therefore have that the integral I 3 is approximately given by 


h{p,k, qi,q2,zi- zq) 


— 2 TTi 
2 tt 


1 (dvra^) 1 /j(u)(f+i:}^)(zi-zo) _ i(u}i+izio){zi-zo)\ 

(wi + Cbm) M 2 M 1 


(dvra. 


'|2 p-U2(z2-zi) 


M2-P+ Ml - M2 


'|2 p-Ul(z2-zi) 




J(cdO+‘^m)(zi-zo) _ p-U2(zi-zo 


+ 


(dvra. 


Ml - M2 


(e 


i{ojO+Cjm)izi-zo) _ -fll{zi-ZO 


(B.d.21) 


This integral reduces to the short path length case and the dimensions are consistent. 
Now, in the well separated case, the last two lines are no longer exponentially suppressed. 
However, in the contact limit, they are equal but opposite in sign and therefore also 
cancel each other. Therefore, for the contact case in which we are interested, they do not 
contribute and we may leave them out of the remainder of the calculation. It is important 
to notice that there is no factor of 1/2. The well separated case is not relevant for the 
short path length limit. 


We therefore have the Matrix element, unchanged in the contact limit 

d^qi 




(27r; 


^.e-*qi-bi^;(0,q,)i;(0,q2)x 


X 

X 


d q 2 -iq,-b 2 n„- 
(271)2® 




e- (k-qi) 


^i(Luo+‘^m)(zi-zo) _ ^i(u}o-cvi)(zi-zo) 


(k — qi)2 + M2a;2 + 

a2[c,ai](ra^ra2) 


(B.d.22) 
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B.4.2 ^ 2 , 0,2 

The result for this diagram is calculated in precisely the manner of ^ 2 , 0,1 with a simple 
change of variable name that involves interchanging all I’s and 2's. Since scattering centres 
are identical then, we must symmetrise these two diagrams which is done by effectively 
multiplying each diagram by 1/2. This result still holds in the short path length limit 
since the change of variables applies across the board. 


^ 2 A 2 j 


d^qi 

( 2 vr )2 


e-'‘ii-^iu(0,qi)u(0,q2)x 


X 

X 


f ^ g-»q2'b22^q 

J ( 27 r) 2 ® 

^i{LUO+id 

771 ){zi-zo) _ g 


e- (k-qi) 


k — q^)2 + 


LJO-OJl){zi-Zo) 


ai [c, aol(r„,,r„,) 


(B.4.23) 


It is also worth mentioning that the diagrams in in Appendix F of [10] will also be 
suppressed in the same manner in the short path length limit. 








Appendix C 


Detailed calculation of Energy 
Loss Formula 


In this section I outline the process whereby one obtains an expression for the first order 
(in opacity) energy loss by using equation (2.3.14). 

We will start with equation 3.0.9 

where is the sum of all diagrams with one interaction with a scattering centre and ^2 
the sum of all diagrams with two interactions with a single scattering centre. Therefore, 
in order to calculate this quantity, we must sum the relevant diagrams first. Once the 
summation is done, we square the amplitude (and take the real part of the trace), average 
over initial and sum over final states. This means that there is an averaging over impact 
parameter (introducing a factor of 1 /^_l) and a sum over scattering centres (introducing 
a factor of A). 

C.l Calculation of ^Tr(|^ip) 

Using the results from sections B.1.1, B.1.3 and B.1.2, given in equations (B.1.8), (B.1.25) 
and (B.1.15) and performing a rearrangement of terms so as to group like phases together. 
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we obtain the sum of the one scattering centre diagrams. 


.^1 + -^ 1 , 1,0 + -^ 1 , 0,1 


=J{p)e^^P+'^>^{-i){2ig,)T, 




v( 0 ,qi)e 


-iqi'bi 




+ 


e • k 

m 2 4 

e • k 

k 2 + 

e • (k - qi^ 


m 2 + k^ + x 2 M 2 L 


“'i (27r)2 
[e*(‘^o+<:'m) _ i]a^c 




COl 


(k — qi )2 + M 2 x 2 + m 2 




(a;o+lI'm) (21-20) _ g*(‘^0-'^l)(2l-2o) 


[c,ai]Ta 


(C.1.1) 


The term that arises from the short path length generalization is the second term on the 
third to last line. It is necessary to reorganise the big bracket - we do this so as to group 
terms with like exponential components. One has to be careful with the color matrices a 
and c as they do not commute 


=-#1,0,0 + -#1,1,0 + -#1,0,1 
=J{p)e^^P+’^><^{-i){2ig,)Ta 


d 2 qi 

(271)2 


i;( 0 ,qi)e 


-iqi'bi 


( e-(k-qi) 

e • k \ 

V (k — q ;^)2 + M 2 x 2 + m 2 

m 2 + k^ + x 2 m 2 / 


J{u)o+Cj 

m )(21-2o) 


c,ai] 


e- (k-qi) 


(k — qi )2 + M 2 x 2 + m 2 

e • k 

aic 


„i(a;o-aJi)(2i-2o) 


C, ai] 


m2 + k^ + x2M2 


1 

+ X 


€ • k 


2 m 2 + k^ + x 2 M 2 


g-Ml(2l-2o)g^^ 


(C.1.2) 
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To simplify the process above, consider the following short hand: 


fk = 


fa — 


e • k 




e- (k-qi) 


(k — q]^)2 + + 777,2 

Wom = (wo - OJm)izi “ ^o) 

Wol = (wo - Wl)(2;i - Zo) 

-UJirn = Wol - Wom = (wq - Wi)(zi - Zq) - (wq - OJm){zi - Zq) 

= -(Wi- 0Jm){zi - ^o) 

as = 

2 


(C.1.3) 


Then we may write 


■/#! —.^ 1 , 0,0 + -^ 1 , 1,0 + >^ 1 , 0,1 

=J{p)e^(P+’^^-°{-i){2igs)Ta, J ^v{0, x 

X {fq - /fc(l - a 3 ))e*‘^°’"[c,a] - [c, a] - fk{l - 03 ) 00 , 


(C.1.4) 


which is a simpler matter to square. The color matrices have to be handled with care; we 
will need a few results: 

[c,a][c,a]* = (-l)[c,a][c,a] 

= — (co — ac){ca — ac) 

= —{caca — caac — acca + acac) 

= 2ccaa — 2caca, (Section A. 8 ) (C.1.5) 


[c, a]ca = (co — ac)ca 
= caca — acca 

= acac — acca, (change of variables in the first term) 

= ac{ac — ca) = ac[a, c] ^ ca[a, c] = [c, a\ac (C.1.6) 


[c, a]CO = {ca — ac)ca = caca — acca 
= —{(P'a^ — caca) 

,(? (xl 


(C.1.7) 


acca = aacc 


(C.1.8) 
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Lastly, we will use the following parameter but our result is slightly more general than 
the result used in [10]; On line (C.1.9), the identity color matrix is assumed in [10] to 
be the identity matrix in the adjoint representation - effectively assuming that the hard 
parton is a gluon. We will remain more general to allow for a later comparison of gluon 
vs. quark energy loss, at the expense of a slightly more complicated expression. 

a = Tr (c^a^ — caca) 

= Tr (c^a^ — caca — if°'^’^cba) 

= Tr(-i/“^"c6a) 

= Tr(-i 1/2 iC 2 {G)aa) 

= ^C2(G)Tr(C2(r)i) (C.1.9) 

= ^C'2(G')C2(r)d(r) (C.1.10) 

Armed with the above results and shorthand, we have compute the sum of the single 
interaction diagrams. At this point we keep, for argument’s sake, the phases. At the end 
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we will examine the result in the limit that //j w*. Noting that 


d'j' 


TV(| = N\J{p)\\Ag, 


2,^ 2^ 1 f d'qi C2iT) 


At J ( 27 r )2 dA 


X |2a(^/g - /fe(l - as)^ (^fq - /fe(l - + 2a/, 


— 2 a 


fq - /fc(l - a3)^/ge*“l"* + (jq - /fc(l - a*^)^ fqe~ 


—a 


{fq - fk{l - a3))/fc(l - 


+ {fq - fk{l - al))fk{l - 
+a(^fqfke^°^{l - ale-^°^) + fqfke-^^°^{l - 


+ Tic^aykfq{l - a3e*‘"°-)(l - 


2,, 2^ 1 f C2{T) 


x< 4 a/^(l -COSWOm) + 2 a/|(l -COSWOm) 


- 4a/fc/g(l - cos Uim) + afqfk2 cos wqi + Tr c^a^f^+ 


+e 


—fiAz 


fl Tr (cos UQm - 1) - Tr c^a^fl cos wom 


2„2 ^2 


+ fkfqOi{cOSLOOm “ COSWoi) 


+ ^/|Tr c^a^e 


C.2 Calculation of Tr(^Q*^ 2 ) 

The averaging over impact parameter plays a crucial simplifying role here, 

4l± 


dij 


(C.1.11) 


(C.1.12) 


(C.2.1) 


So, consider the results in equations (B.2.30), (B.3.29), (B.3.62), (B.4.22) and (B.4.22). 
If they are added up and rearranged slightly, again with the purpose of grouping terms 
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with similar exponential components, one finds that 


2 , 0,3 + ^^ 2 , 0,0 + -^ 2 , 2,0 + ^ 2 , 0,1 + -^ 2 , 0,2 


k 2 + M 2 x 2 + m 2 


i{LLJCl+U) 

m )(2l-2o) 


+ 


e • k 


+ M 2 x 2 + m 2 

e • k 

k^ + M2 x 2 + m2 


a2aic(e*("'°-“(i2))(^i-^o) _ i) 


casai ei(‘^o-^(i2))izi-zo) 


_l_ g-m(^i-^o) _ Ig-m(^i-zo)^^ 

€ • (k - q^) , 

^(k-qi)2 + M2x2 + m2“2[^’ 


i(u)(i+Cj 

m )(21-2o) 


_ gi(a;o-wi)( 2 i- 2 o)^ 


— e 


i(a;o-wi)( 2 i- 2 o 




(C.2.2) 


Here we have used the fact that 

2 2,2 2,2 2 

Ati = qi + M = q 2 + /^ = /^l- 


Now we need , which is given in [10] to be 


^0 ~ 


(-2i5s) 


e • k 

k^ + m2 + M2 x2 


giwo^Og 


Using the same shorthand as before, we may write 
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{^^^ 2 )= N\J{p)\\Agl) — 


d\i C2{T) 


At J (27r)2 dA 

/I (2a cos ojom — 2a — Tr c^a^) 
+ 2afkfq{cOSUJ0m “ COSWqi) 

+ e~^^''fk Tr cos wom - 


(C.2.3) 


C.3 Addition of summed terms 


It now remains simply to add the contributions from equations (C.1.12) and (C.2.3) which, 
after a number of cancellations, gives 

^^'(1 

X - 2fq{fk - fq){l - COSWlm) 

+e“^^^^a/|(l - coswom) + /|(coswom - 1) Tr 

+/fc/qa(coswom - COSWoi) ) (C.3.1) 


C.4 Final Calculation 

We now want to extract dEj^^^ = ud^Ng. From [10], we know that 


which we can use, along with equation (3.0.9), (equation (10) from [10] to find 

= Aj- + |-»TV(.^„V^2>) (C-4.2) 

We now do a kind of averaging of the position of the zi because the fraction of energy 
carried off by the radiated gluon seems to depend on the distance between the production 
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{zq) and the interaction with the medium (zi): 

/ dE\ \ dE f dE f (zi — zq) \ 2 

This is an assumed distribution. It tries to take into account the rapidly exopanding 
medium so that the distance between scattering centres increases exponentially. Noting 
that 


= dk^d^lL = 

= d{xE~^)d^'k = E'^dxd^'k = 2Edxd^\L 


For clarity, let 


★ = 


-‘^fqifk - /g)(l - COSWlm) 


+e - coswom) + /|(cosa;om - 1) Tr 

+ fkfqa{cOS0J0m “ COSWqi) ) 


(C.4.4) 


We will also be using a ‘clever “1” From [12], we have a number of relations that will 
be useful. 

• For a homogeneous rectangular target of thickness L, the density is p = N/LA±. 

• Since the mean free path is A = l/po-^i, the opacity is simply L/X = Naei/A±. 

• In perturbation theory, p/A Rs dvra^p. 

And so we can wrangle a convenient ‘one’: 

— = 47ra„p 

A 

Aira^N 
~ LA± 

^ ^ p‘^LA± 

XAira'^N 


(C.4.5) 
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We start with 








X 


(27r)3 2 

xE+/ J dzip{zi)[d.] 

N^i^nas) 1 /• 2 „ (47ra,)2 

xE+y J dzip{zi) 

Here we multiply by the above strategic ‘1’ and divide dx out 


(C.4.6) 


dEll\_ME+Cn f. (47ra,)3 


dx 


d qi- 


X / d^k 


, 2 ,^ 1 p^L^ 1 


(27r)^ XAtto^M (27r)3 


dzip{zi) H* 


(C.4.7) 


= / d / d k 


, 2 , -B“''C'ij647r3Q;3^2^ 


(qf + /u2)2327r®A4a2 


X J dzi — exp 


2htp>K*i 


Using the fact that E~^ ~ 2E 


CrUsLE f d^qi pf 


d^k 


ttA 


TT {p? + qf47r 


dzi exp 


-2Up>}|M. 


which expands to 


(C.4.8) 


CRagLE f d‘^c{i p^ d^k 
■ X 


ttA 7 tt {p^ + qfTT 

- 2fq{fk - fg){l - COSWlm) 






+ xe 


—fiAz I f2 


/|(1 -COSWOm)^ - 


Tr 


a 


+ fkfqOt{ cos WOm “ COS Wol) 


(C.4.9) 
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At this point we pause to look at a few results. Firstly, the zi integral for the first term 
in the square brackets (the leading order term derived by [10] is analytically known: 


dzi exp \ —2 


, (^1 - ^o) 1 2 


L j L 


y(l - COS {(CJI + Cjm){zi - 20 )}) 


(cJl + U}rn)‘^L‘^ 

4 + (oJi + 

— qi)^ + ^ 

^ [(k - qi)2 + m2 + M2 x2]^ + 


(C.4.10) 


The zi integral can also be performed analytically for the next to leading order terms. 
They are cumbersome, but their form is insightful so I present them here. Term by term 
then, we have, for the next to leading order correction terms for the short separation 
distance generalization. 


^ dziexp I - 2 ^^^ ^ |^(cos{(cuo - (Lm){zi - zq)} - 1) 

2 / 


2L (tCg d)rn) 


(2 + L/ri) (4 + 4L/ii + L‘^{{loq — cDm)^ + A^i)) 

X (cos{(a;o - OJmJizi - Zq)} - cos{(wo - (Zli){zi - Zo)}) 
2(2 + Ljj,\) 


(C.4.11) 


+ 


4 + 4L/ii + LP‘{{ijJq - wi)2 + gi\) 

2(2 + T/ii) 

4 + ALni + L'^ifujQ — CbmY' T Li) 


(C.4.12) 


Although further simplifications can be performed on the leading order terms and despite 
the fact that analytical expressions exist for the integrals over the separation distance, 
the expression becomes unnecessarily involved. I therefore present here the full formula 
with corrections, but before the evaluation of the zi integral: 
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^ ^( 1 ) CnasLE f 




dAzp{Az) X 


J TT (^2 _|_ q2^2 4^ 
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